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Abstract 

We perform complete group classification of the general class of quasi linear wave equations 
in two variables. This class may be seen as a broad generalization of the nonlinear d'Alembert, 
Liouville, sin/sinh-Gordon and Tzitzeica equations. In this way we derived a number of new 
genuinely nonlinear invariant models with high symmetry properties. In particular, we obtain four 
classes of nonlinear wave equations admitting five-dimensional invariance groups. Applying the 
symmetry reduction technique we construct multi-parameter families of exact solutions of those 
wave equations. 



Introduction 
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Cp ■ More than century ago Sophus Lie introduced the concept of continuous transformation group into 

mathematical physics and mechanics. His initial motivation was to develop a theory of integration 
of ordinary differential equations enabling to answer the basic questions, like, why some equations 
are integrable and others are not. His fundamental results obtained on this way, can be seen as a 
. far reaching generalization of the Galua's and Abel's theory of solubility of algebraic equations by 

^ I radicals. Since that time the Lie's theory of continuous transformation groups has become applicable 

to an astonishingly wide range of mathematical and physical problems. 

It was Lie who was the first to utilize group properties of differential equations for constructing of 
I their exact solutions. In particular, he computed the maximal invariance group of the one-dimensional 

heat conductivity equation and applied this symmetry to construct its explicit solutions. Saying it 
the modern way, he performed symmetry reduction of the heat equation. Since late 1970s symmetry 
reduction becomes one of the most popular tools for solving nonlinear partial differential equations 
(PDEs). 

By now symmetry properties of the majority of fundamental equations of mathematical and the- 
oretical physics are well known. It turns out that for the most part these equations admit wide 
symmetry groups. Especially this is the case for linear PDEs and it is this rich symmetry that enables 
developing a variety of efficient methods for mathematical analysis of linear differential equations. 
However, linear equations give mathematical description of physical, chemical or biological processes 
in a first approximation only. To provide a more detailed and precise description a mathematical 
model has to incorporate nonlinear terms. Note that some important equations (for example, the 
Yang-Mills equations) of theoretical physics are essentially nonlinear in a sense that they have no 
linearized version. 
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Hyperbolic type second-order nonlinear PDEs in two independent variables play a fundamental 
role in modern mathematical physics. Equations of this type are utilized to describe different types of 
wave propagation. They are used in differential geometry, in various fields of hydro- and gas dynamics, 
chemical technology, super conductivity, crystal dislocation to mention only a few applications areas. 
Surprisingly the list of equations utilized is rather narrow. In fact, it is comprised by the Liouville, 
sine/sinh-Gordon, Goursat, d'Alembert and Tzitzeica equations and a couple of others. Popularity of 
these very models has a natural group-theoretical interpretation, namely, all of them have non-trivial 
Lie or Lie-Backlund symmetry. By his very reason some of them are integrable by the inverse problem 
methods (see, e.g.,[2-|ni) or are linearizable [3~|S1) and completely integrable |3IH]- 

In this connection it seems a very important problem to select from the reasonably extensive class 
of nonlinear hyperbolic type PDEs those enjoying the best symmetry properties. Saying 'reasonably 
extensive" we mean that this class should contain the above enumerated equations as particular cases 
on the one hand, and on the other it should contain a wide variety of new invariant models of potential 
interest for applications. The list of the so obtained invariant equations will contain candidates for 
realistic nonlinear mathematical models of the physical and chemical processes enumerated above. 

The history of group classification methods goes back to Lie itself. Probably, the very first paper 
on this subject is P], where Lie proves that a linear two-dimensional second-order PDE may admit at 
most a three-parameter invariance group (apart from the trivial infinite-parameter symmetry group, 
which is due to linearity). 

The modern formulation of the problem of group classification of PDEs was suggested by Ovsyan- 
nikov in ^U]. He developed a regular method (we will refer to it as to Lie-Ovsyannikov method) for 
classifying differential equations with non-trivial symmetry and performed complete group classifica- 
tion of the nonlinear heat conductivity equation. In a number of subsequent publications more general 
types of nonlinear heat equations were classified (review of these results can be found in ) . 

However, even a very quick analysis of the papers on group classification of PDEs reveals that 
the overwhelming majority of them deal with equations whose arbitrary elements (functions) depend 
on one variable only. There is a deep reason for this fact, which is that Lie-Ovsyannikov method 
becomes inefficient for PDEs containing arbitrary functions of several variables. To achieve a complete 
classification one either needs to specify the transformation group realization or somehow restrict 
arbitrariness of functions contained in the equation under study. 

Recently, we developed an efficient approach enabling to overcome this difficulty for low dimensional 
PDEs jl2l I13j . Utilizing it approach we have obtained a final solution of the problem of group 
classification of the general quasi-linear heat conductivity equations in two independent variables. 

In this paper we apply the approach mentioned above to perform group classification of the most 
general quasi-linear hyperbolic type PDE in two independent variables. 

1 Group classification algorithm 

We begin this section by formulating the problem to be solved. Then we present a brief review of the 
already known results. Finally we give a short description of our approach to group classification of 
PDEs (for the detailed account of the necessary facts, see (^31))- 

While classifying a given class of differential equations into subclasses, one can use different clas- 
sifying features, like linearity, order, the number of independent or dependent variables, etc. In group 
analysis of differential equations the principal classifying features are symmetry properties of equa- 
tions under study. This means that classification objects are equations together with their symmetry 
groups. This point of view is based on the well-known fact that any PDE admits a (possibly trivial) Lie 
transformation group. And what is more, any transformation group corresponds to a class of PDEs, 
which are invariant under this group. So to perform group classification of a class of PDEs means de- 
scribing all possible (inequivalent) pairs (PDE, maximal invariance group), where PDE should belong 
to the class of equations under consideration. 



We perform group classification of the following class of quasi-linear wave equations 



Utt = Uxx + F{t,X,U,Ux). (1-1) 



Here F is an arbitrary smooth function, u = u{t,x). Hereafter we adopt notations ut = Ux = 

, Utt - Q^, 

Our aim is describing all equations of the form Hl.l|) that admit nontrivial symmetry groups. The 
challenge of this task is in the word all. If, for example, we fix the desired invariance group, then the 
classification problem becomes a University textbook exercise on Lie group analysis. A slightly more 
cumbersome (but still tractable with the standard Lie-Ovsyannikov approach) is the problem of group 
classification of equation with arbitrary functions of, at most, one variable. 

As equations invariant under similar Lie groups are identical within the group-theoretic framework, 
it makes sense to consider non similar transformation groups |14| I15j only. The important example of 
similar Lie groups is provided by Lie transformation groups obtained one from another by a suitable 
change of variables. Consequently, equations obtained one from another by a change of variables have 
similar symmetry groups and cannot be distinguished within the group-theoretical viewpoint. That is 
why, we perform group classification of (|l.lj) within a change of variables preserving the class of PDEs 

CU). 

The problem of group classification of linear hyperbolic type equation 

utx + A{t, x)ut + B{t, x)ux + C{t, x)u = (1.2) 

with u = u{t,x), was solved by Lie [HI (see, also, [El). In view of this fact, we consider only those 
equations of the form Hl.l|) which are not (locally) equivalent to the linear equation (|1.2|) . 

As we have already mentioned in Introduction, the Lie-Ovsyannikov method of group classification 
of differential equations has been suggested in jlUj . Utilizing this method enabled solving group 
classification problem for a number of important one-dimensional nonlinear wave equations: 



Utt = Uxx + F{u) 
Utt = [f{u)Ux]x 
Utt = f{Ux)Uxx 
Utt = F{Ux)Uxx + H{Ux) 
Utt = F{Uxx) 
Utt = u'^Uxx + f{u) 
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1201-122] 
|22l 123 

m 
m 
m 



Utt + f{u)ut = {g{u)ux)x + h{u)ux. Uni 

Analysis of the above list shows immediately that arbitrary elements (= arbitrary functions) depend 
on one variable, at most. This is not coincidental, indeed, the Lie-Ovsyannikov approach works 
smoothly for the case when the arbitrary elements are functions of one variable. The reason for this 
is that the obtained system of determining equations is still over-determined and can be effectively 
solved within the same ideology used while computing maximal symmetry group of PDEs containing 
no arbitrary elements. 

The situation becomes much more complicated for the case when arbitrary elements are functions 
of two (or more) arguments. By this very reason group classification of nonlinear wave equations 



Utt = {f{x,u)Ux)x] W\ 

Utt + >^Uxx = g{u, Ux); UHl HH] 

Utt = [f{u)ux + g{x, u)]x; 1301 - 

Utt = f{x,Ux)uxx + gix,Ux) 1331 



is not complete. 

We suggest new efficient approach to the problem of group classification of low dimensional PDEs 
in |12|ll,Sj. This approach is based on Lie-Ovsyannikov infinitesimal method and classification results 
for abstract finite-dimensional Lie algebras. It enables us to obtain the complete solution of group 
classification problem for the general heat equation with a nonlinear source 



Later on, we perform complete group classification of the most general quasi-linear evolution equation 



We use the above approach to achieve a complete solution of the group classification problem for the 
class of equations 

Our approach to group classification of the class of PDEs consists of the following steps (for 
more details, see [36]): 

L Using the infinitesimal Lie method we derive the system of determining equations for coefficients 
of the first-order operator that generates symmetry group of equation (Note that the 

determining equations which explicitly depend on the function F and its derivatives are called 
classifying equations). Integrating equations that do not depend on F we obtain the form of the 
most general infinitesimal operator admitted by equation Hl.lj) under arbitrary F. Another task 
of this step is calculating the equivalence group £■ of the class of PDEs (|1.H) . 

II. We construct all realizations of Lie algebras of the dimension n < 3 in the class of operators 
obtained at the first step within the equivalence relation defined by transformations from the 
equivalence group £. Inserting the so obtained operators into classifying equations we select 
those realizations that can be symmetry algebras of a differential equation of the form (|1.H) . 

III. We perform extension of the realizations constructed at the previous step to realizations of 
higher dimensional (n > 3) Lie algebras. Since extending symmetry algebras results in reducing 
arbitrariness of the function F, at some point this function will contain either arbitrary functions 
of at most one variable or arbitrary constants. At this point, we apply the standard classification 
method, which is due to Lie and Ovsyannikov, to derive the maximal symmetry group of the 
equation under study thus completing its group classification. 

As a result of performing the above enumerated steps we get the complete list of inequivalent 
equations of the form 1)1.11) together with their maximal (in Lie's sense) symmetry algebras. 
We say that the group classification problem is completely solved when it is proved that 

1 ) The constructed symmetry algebras are maximal invariance algebras of the equations under 

consideration; 

2 ) The list of invariant equations contains only inequivalent ones, namely, no equation can be 

transformed into another one from the list by a transformation from the equivalence group £. 

2 Preliminary group classification of equation (1.1) 

According to the above algorithm we are looking for infinitesimal operator of symmetry group of 
equation ()1.1() in the form 



Ut = Uxx + F{t,X,U,Ux). 



m-m 



Ut = f{t, X, U, Ux)u^x + git, X, U, Ux). 



Q = r(t, X, u)dt + i{t, X, u)dx + r/(t, x, ■u)5, 




where r, ^, rj are smooth functions defined on an open domain ft of the space 1/ = x of indepen- 
dent M? = (t, x) and dependent = (n) = u(t, x) variables. 

Operator (|2.1|) generates one-parameter invariance group of iff its coefficients r, ^, ry, e satisfy 
the equation (Lie's invariance criterion) 



^tt _ ^xx _ _ _ ^p^ _ ^xp^^ 



= 0, (2.2) 

o 



where 



yp* = Dt{v)-UtDt{T)-uMO, 

if"" = D^{7]) - utD^ir) - u^D^{C), 

if'' = Dt{^')-uuDt{T)-Ut,Dt{0, 

ip^^ = D^iip"") - UtxD^ij) - u^xDa,(0 

and Dt,Dx are operators of total differentiation with respect to the variables t,x. As customary, 

by writing we mean that one needs to replace uu and its differential consequences with the 

expression Uxx + F and its differential consequences in H2.2() . 

After a simple algebra we represent ()2.2I) in the form of system of four PDEs: 

(1) ^„ = T„ = rjuu = 0, 

(2) Tt-^x=0, ^t-Tx = 0, 

(3) 27]tu + TxFu, = 0, (2.3) 

(4) r]tt - r]xx - 'iUx'Tlxu + [riu - 2rt]F - rFt - ^F^ 
-rjFu - [-qx + Uxitlu - Cx)]Fux = 0. 

It follows immediately from (1) that r = T{t,x), ^ = ^{t,x),rj = h{t,x)u + r{t,x). In the sequel 
we differentiate between the cases Fu^u^ / and Fu^u^ = 0- 

Case Fu^ux 0- It follows from (3) that = ht = 0. Taking into account this fact and equation 
(2) we obtain r = Ai + Ai, = Ax + A2, h = h{x), where A, Ai, A2 are arbitrary real constants. 

Case Fu^ux = 0. If this is the case, then F = g{t,x,u)Ux + f{t,x,u), where / and g are arbitrary 
smooth functions. 

Given the condition 7^ 0, it follows from (3) that = ht = 0. So that taking into account 
equation (2) we arrive at the expressions r, ^, h given in the previous case. 
If gu = 0, then 7^ (otherwise equation (|1.H1 becomes linear). 
Let g = 0. It follows from (3) that rj = h{x)u + r{t,x). Equation (4) now reads as 

nt - Txx - 2h'ux + [h- 2Tt]f - Tft - ifx - [hu + r]fu = 0. 

As functions r, ^, h, r, f do not depend on Ux, h' = 0. Hence r] = mu + r{t, x), where m is an arbitrary 
real constant. 

Furthermore, if 5 = g{t, x) ^ 0, then it is straightforward to verify that system of equations ((3), 
(4)) is equivalent to the following equations: 

2ht = -Txg, 2hx = -ng - rgt - Cgx, 

{hu - hxx)u + rtt - rxx + f[h - 2Tt] - rft - Cfx - [hu + r]/„ - {hxU + rx)g = 0. 

Integrating (2) yields r = ip{9) + ^ = —^{0) + ipic), where ipjip are arbitrary smooth functions 

oi6 = t — x,a = t + x and we arrive at the following theorem. 



Theorem 1 Provided F^^ux 7^ 0; the maximal invariance group of equation il.lp is generated by the 
following infinitesimal operator: 

Q = {Xt + Xi)dt + {Xx + X2)d^ + [h{x)u + r{t, x)]du, (2.4) 

where A,Ai,A2 are real constants and h = h{x), r = r{t,x), F = F{t,x,u,Ux) are functions obeying 
the constraint 

d^h dh 
nt - Txx - -r^u - 2—Ux + {h- 2X) F 
ax^ dx 

-(At + Ai) Ft - {Xx + Aa) F^ - {hu + r) F^ (2.5) 
-{tx + + {h- X)ux) Fu^ = 0. 

If F = g(t, x,u)ux + f{t,x,u), gu ^ 0, then the maximal invariance group of equation 
is generated by infinitesimal operator i2.4\ l, where A,Ai, A2 are real constants h,r,g,f are functions 
satisfying system of two equations 

-2h' - Xg = (At + Xi)gt + (Ax + X2)gx + {hu + r)gu, 

-h"u + rtt - Txx + ih- 2X)f = (At + Ai)/t + (Ax + A2)/. (2.6) 
+ {hu + r)fu + g{h'u + rx). 

Next, if F = g{t,x)ux + f{t,x,u), q ^ 0, f^u 7^ 0, then the infinitesimal operator of the invariance 
group of equation hl.l\) reads as 

Q = T{t, x)dt + i{t, x)dx + {h{t, x)u + r{t, x))du, 

where r, ^, h, r, g, f are functions satisfying system of PDEs 

Tt-L = 0, it-Tx = 0, 
2ht = -Txg, 2hx = -ng - rgt - Cfe, 
{hu - hxx)u + rtt - Txx + f{h - 2Tt) - rft 
-ifx - {hu + r)fu - {hxU + rx)g = 0. 

Finally, if F = f{t,x,u), fuu 7^ 0, then the maximal invariance group of equation is generated 

by infinitesimal operator 

Q = [ip{e) + ^P{a)]dt - [ip{e) - i;{a)]dx + [ku + r{t, x)]du, 
where /cGM, 6 = t — x, a = t + x and functions tp, ip, r, f and constant k satisfy the following equation: 

rtt - Txx + [k- 2ip' - 2V']/ -((/? + ^)/t + 
+(</'- ^)/x-(A;n + r)/„ = 0, </'' = ^, ^' = ^- 

Summing up the above considerations we conclude that the problem of group classification of 
equation reduces to the one of classifying equations of more specific forms 

Utt = Uxx + F{t , X, 1/, tlx 

), F^.u.y^O; (2.7) 

utt = Uxx + 9{t,x,u)ux + f{t,x,u), 9« / 0; (2.8) 

utt = Uxx + 9{t,x)ux + f{t,x,u), 5 / 0, /un / 0; (2.9) 

Utt = Uxx + f{t,x,u), fuu^O- (2.10) 



Consider the last two equations. By the change of variables 



t = t — x, x = t + x, u = v{t,x) 

we reduce them to equations 

Vtx = ^f{t,x,v), 

Vtx = -^9{t,x){vt-vsi) + ^f{t,x,v). (2.11) 
Now, making the change of variables 

i = t, X = x, v{i, x) = A{t, x)v, 
where A = exp [— ^ J g{t, x)dx] , we transform ()2.11|1 to become 

Vix = Q<? - A"^Ai^ - ^gA-^AiV + ^A-^A^gv + A'V. 
Hence we conclude that the following assertion holds true. 

Assertion 1 The problem of group classification of equations h2.!J\) . h2.1(l\) is equivalent to the one of 
classifying equations 

utx = git,x)ux + fit,x,u), gx ^ 0, fuu ^ 0; (2.12) 
utx = fit,x,u), fuu^O. (2.13) 

Note that condition gx ^ is essential, since otherwise (|2.12|) is locally equivalent (|2.13|) . 

Summing up, we conclude that the problem of group classification of ()1.1|) reduces to classifying 
more particular classes of PDEs (|2.7)) . (|2.81 (|2.12)) . (|2.13|) . In what follows, we provide full calculation 
details for equations (|2.9|) and 1)2.10(1 only. The reason is just to save space and still be able to present 
all details of the algorithm. 

First, we consider equations (|^ . (|TT^ . (|Tn|) . 



3 Group classification of equation (2.8) 

According to Theorem ^ invar iance group of equation (|2.8() is generated by infinitesimal operator (|2.4j) . 
And what is more, the real constants A, Ai,A2 and functions h,r,g,f satisfy equations H2.6|) . System 
(|2.6|) is to be used to specify both the form of functions /, g from (|2.8|) and functions h, r and constants 
A, Ai,A2 in ()2.4|) . It is called the determining (sometimes classifying) equations. 

Efficiency of the Lie method for calculation of maximal invariance group of PDE is essentially 
based on the fact that routinely this system is over-determined. This is clearly not the case, since we 
have only one equation for four (!) arbitrary functions and three of the latter depend on two variables. 
By this very reason the direct application of Lie approach in the Ovsyannikov's spirit is no longer 
efficient when we attempt classifying PDEs with arbitrary functions of several variables. 

Next, we compute the equivalence group £ of equation (|2.8j) . This group is generated by invertible 
transformations of the space V preserving the differential structure of equation (|2.8() (see, e.g., [Tlj). 
Saying it another way, group transformation from £ 

D(i,x,v) 

t = a{t,x,u), x = (3{t,x,u), v = U{t,x,u), ' — ^ — r 7^ 0, 

J-^ {l J X J 111 ] 



should reduce (|2.8jl to equation of the same form 

vu = Vxx + g{i, X, v)vx + fit, x,v), gv 7^0 

with possibly different f,g. 

As proved by Ovsyannikov |14J, it is possible to modify the Lie's infinitesimal approach to calculate 
equivalence group in essentially same way as invariance group. We omit the simple intermediate 
calculations and present the final result. 

Assertion 2 The maximal equivalence group E of equation 112.^) reads as 

t = kt + ki, x = ekx + k2, v = X{x)u + Y{t, x), (3-1) 
where A; / 0, X / 0, e = ±1, k, ki,k2 € and X, Y are arbitrary smooth functions. 
This completes the first step of the algorithm. 

3.1 Preliminary group classification of equation (2.8). 

First, select equations of the form (|2.8jl admitting one-parameter invariance groups. 

Lemma 1 There exist transformations that reduce operator \2.4\) to one of the six forms: 

Q = m{tdt + xdx), m 7^ 0; Q = dt + pd^, /3 > 0; 

Q = dt + a{x)udu, (J ^ 0; Q = d^; (3.2) 
Q = a{x)udu, Q = e{t,x)du, 9^0. 

Proof. Change of variables reduces operator (|2.4j) to the form 

Q = k{\t + \i)di + ek{\x + \2)dx + [Yt{\t + Ai) + {\x + \2){X'u + Y^) + X{hu + r)]d.,. (3.3) 

If A 7^ in (|2.4() . then putting ki = X^'^Xik, k2 = eA^^A2A;, and taking as X,Y {X ^ 0) integrals of 
system of PDEs 

X'{\x + \2) + Xh = 0, 

Yt{\t + Ai) + Yx{\x + As) + Xr = 

in (|3.1|1 . we reduce (|3.3|1 to the form 

Q = X{tdt + xds). 
Provided A = and Ai / 0, we analogously obtain 

Q = dt + f3ds, /3>0; Q = dt + a{x)vd^, a / 0. 

Next, if A = Ai = 0, A2 / 0, in (|^ . then putting k = eX^^, and taking as X, Y (X 7^ 0) 
integrals of equations 

X2X' + hX = 0, Yx + rX = 0, 

we reduce operator ()3.3|) to become Q = dx. 

Finally, the case A = Ai = A2 = 0, gives rise to operators Q = a{x)vdv, Q = 6{t,x)dy. 
Rewriting the above operators in the initial variables t, x completes the proof. 



Theorem 2 There are exactly five inequivalent equations of the form that admit one-parameter 

transformation groups. They are listed below together with one- dimensional Lie algebras generating 
their invariance groups (note that we do not present the full form of invariant PDEs just the functions 
f and g) 

A\ = {tdt + xdx) ■■ g = x~'^g{i^, u), 

f = x~'^f{i^, u),'il^ = tx-^,gu / 0; 

Al = {dt + Pd^) : g = g{r],u), f = f{r],u), 

7j = x-pt, P>0, gu^O; 
= {dt + (T{x)udu) : g = -2a' ln|n| +g{p,x), 

f = {a' a^^Yulv? \u\ — a'a^^g{p, x)u In \u\ — a^^a"u In \u\ + uf{p, x), 

/9 = nexp(— tfj), cr 7^ 0; 
^1 = (dx) ■■ g = g{t,u), f = f{t,u), gu + 0; 

A\ = {a{x)udu) : g = -2a' cf-^ In \u\ + g{t, x), f = {a'a'^ fu In^ |n| 
-{a-^a" + a-^ a' g{t,x))u In \u\ + uf{t,x), a' / 0. 

Proof. If equation H2.8|) admits a one-parameter invariance group, then it is generated by operator 
of the form (|2.4j) . According to Lemma ^ the latter is equivalent to one of the six operators (|H.2|) . 
That is why, all we need to do is integrating six systems of determining equations corresponding to 
operators (|2.H|) . For the first five operators solutions of determining equations are easily shown to 
have the form given in the statement of theorem. 

We consider in more detail the operator Q = 9{t,u)du. Determining equations (|2.6|) for this 
operator reduce to the form 

Gtt — Oxx = Gfu + 0x9, Ogu = 0, 

whence we get gu = 0. Consequently, the system of determining equations is incompatible and the 
corresponding invariant equation fails to exist. 

Non equivalence of the invariant equations follows from non equivalence of the corresponding 
symmetry operators. 

The theorem is proved. 

Lemma 2 There are no realizations of semi-simple Lie algebras generated by operators of the form 




Proof. To prove the lemma it suffices to check that there are no realizations of the lowest order 
simple Lie algebras by operators (|2.4() . The commutation relations defining those read as j37j : 

so(3) = (61,62,63) : [61,62] =63, [61,63] = -62, [62,63] =61; 

S/(2,M) = (61,62,63) : [6l,62] = 262, [61,63] = -263, [62,63] = 61. 

We start by noting that one of the basis operators 61,62,63 can be reduced to one of the five 
operators (|.'-{.2jl (see. Lemma P). 

We consider in detail the case of operator 

tdt + xd^ (3.4) 
only, since the remaining operators are treated in a similar way. 



Let the basis operator ei of the algebras so(3) and s/(2,M) be of the form (|3.4j) . Computing the 
commutator of ei and Q of the form (|2.4|) yields the relation 

[ei,Q] = -Xidt - X2dx + [xh'u + xr^ + trt]du- 

To satisfy the first two commutation relations for each of the algebras under study, the basis operators 
62, 63 are to be of the form 

aidt + a2dx + {l{x)u + ^(t, x))du, 

where ai,a2 G I^j 7 and /i are smooth functions. It is straightforward to verify that these operators 
cannot satisfy the third commutation relation for either algebra ,3/(2, M) and so(3). 
The lemma is proved. 

Theorem 3 There are no nonlinear equations h2.t^) whose invariance algebras are isomorphic to 
semi-simple Lie algebras or contain those as sub-algebras. 

Proof. Suppose the inverse. Let (|2.8j) be an equation whose invariance algebra contain sub-algebra 
that is semi-simple Lie algebra L. Then by properties of semi-simple Lie algebras, there exist linear 
combinations of the basis elements of L forming the basis of either so(3) or s/(2,M). However, due to 
Lemma 121 there are no realizations of the algebras so(3), s/(2,M) by operators (|2.4)1 . We arrive at the 
contradiction which proves the theorem. 

It follows from Theorem |21 and Levi-Maltsev theorem (see, e.g., |37| I38j ) that nonlinear equations 
1)2. 8|) can admit invariance algebras of the dimension higher than one, provided, (1) those are isomor- 
phic to real solvable Lie algebras, or (2) their finite dimensional sub-algebras are real and solvable. 
Using this fact and also the concept of compositional row for solvable Lie algebras we may perform 
hierarchical classification of invariant equations starting from the lowest dimensional solvable Lie al- 
gebras and increasing dimension by one till we exhaust all possible invariant equations. We start by 
considering two-dimensional solvable Lie algebras. 

There exist two inequivalent two-dimensional solvable Lie algebras |38| 139) 

A2.1 = (61,62) : [61,62] = 0; 
A2.2 = (61,62) : [61,62] = 62. 

To construct all possible realizations of the above algebras we take as the first basis element one of the 
realizations of one-dimensional invariance algebras obtained above. The second operator is looked for 
in the form ()2.4)) . In the case of commutative algebra A2.1 there is no difference between operators 61 
and 62, while for the algebra A2.2 those two operators need separate analysis. We give full computation 
details for the case, when one of the basis elements is of the form A\ given in Theorem [21 

Algebra ^2,1. Let operator 61 be of the form ()3.4() and operator 62 read as 1)2. 4|) . Then it follows 
from the relation [61, 62] = that Ai = A2 = xh' = 0, trt + xr^ = 0. Consequently, we may choose the 
basis elements of the algebra realization in the form {tdt+xdx, {mu-'rr{ip))du) , where m S M, ip = tx~^. 
Provided m = 0, the operator 62 becomes r{'ip)du. As established earlier, this realization does not 
satisfy the determining equations. Hence, m 7^ 0. Making the change of variables 

t = t, x = x, V = u -\- m~^r{'ip) 

reduces the basis operators in question to the form td]; + xdx , mvdy . That is why we can restrict our 
considerations to the realization {tdt + xdx, udu). 

The second determining equation from 1)2.6(1 after being written for the operator udu takes the 
form ugu = 0, whence it follows that this realization does not satisfy the determining equations. So 
the realization A\ cannot be extended to a realization of the two-dimensional algebra ^2.1- 

Algebra ^2.2- If operator 61 is of the form ()3.4)) . then it follows from [61,62] = 62 that A = Ai = 
A2 = 0, xh' = h, trt + xrx = r. 



Next, if 62 reads as (|3.4|) . then we get from [ei, 62] = 62 the erroneous equahty 1 = 0. 

So the only possible case is when 62 = {mxu + xr{ip))du, m 7^ 0, -0 = tx~^, which gives rise to the 
following realization of the algebra j42.2: {tdt + xdx,xudu)- This algebra is indeed invariance algebra of 
an equation from the class H2.8|) and the functions / and g read as: g = —2x~^ In \u\ + x~^g{Tp), f = 
x~^ulii^ \u\ — x~'^g{^l^)uln \u\ + x~^uf{'ifj), -0 = tx~^. 

Analysis of the remaining realizations of one-dimensional Lie algebras yields ten inequivalent A2.1- 
and ^2.2-invariant equations (see the assertions below). What is more, the obtained (two-dimensional) 
algebras are maximal symmetry algebras of the corresponding equations. 

Theorem 4 There are, at most, four inequivalent A2.i-'i'nvariant nonlinear equations Below 
we list the realizations of A2.1 and the corresponding expressions for f and g. 

1) {dt,a{x)udu,) : g = -2a'a^^ln\u\, 

f = {a'a-^fuln^ \u\ -o-~Vnln|u| +uf{x), a' / 0; 

2) {dt, d:r) ■■ g = g{u), f = f{u), Cju + 0; 

3) [d^.dt + n5n) : g = g{uj), f = exp(t)/((j), uj = exp(-t), g^ 7^ 0; 

4) {a{x)udu, dt — -kcr{x)ip{x)udu) ■ g = —2a'a~^ In \u\ + kt + g{x), 



f = {a'a ^)'^uln'^ \u\ — a ^a"uln\u\—a ^a' {kt + g{x))uln\u\ 
+u 



ha'a^H + ^kH^ + ^kg{x) + /(x) 



/c / 0, 0-' / 0, V' = j C7~^dx. 



Theorem 5 There exist, at most, six inequivalent A2,2-in'Variant nonlinear equations \2. Below 
we list the realizations of A2.1 and the corresponding expressions for f and g. 

1) {tdt + xdx, k~^\x\^udu) ■ g = x~^{—2k In |m| -|- ^(V')), 

/ = x~'^u{-k'^ In^ |u| + kg{i{j) In \u\ + k{k - 1) In \u\ + /(V')), 
A; / 0, -0 = te-^ 

2) {dt + ^dx, exp{f3-^x)udu) : g = -2(5-^ In \u\ + g{rj), 

f = p-^uln^ \u\ - {f3-^ + (3-^g{7]))uln \u\ + n/(r/), 
(3 > 0, rj = x- [3t; 

3) {-tdt - xdx, dt + (idx) : g = r^m, f = r^~f{y), /3 > o, 

r] = X- pt, 5„ / 0; 

4) {-tdt - xdx, dt + rnx'^udu) : g = x^'^{2mip + g{u})), 

f = x~^[—2m'ipu — 2mip — 2 — g{uj) + exp{mip)g{uj)], 
m > 0, UJ = uexp{—m7p), ip = tx~^, g^ ^ 0; 

5) {dx,e''udu) : g = -21n|ti| +g{t), f = u\v? \u\- 

-u\n\u\{l + ~g{t)) + uf{t)- 

6) {-tdt - xdx, dx) ■g = t-^g{u), f = t-'^f{u), g^ / 0. 



3.2 Completing group classification of (|2.8|) . 

As the invariant equations obtained in the previous subsection contain arbitrary functions of, at most, 
one variable, we can use the standard Lie-Ovsyannikov approach to complete group classification of 
(|2.8|) . We give the computation details for the case of the first A2.i-invariant equation, the remaining 
cases are treated in a similar way. 

Putting g = —2a'a~^ ln\u\, f = {a'a^^)uln'^ \u\ — a^^a"uln\u\ -\- uf{x), a = cr{x), a' / we 
rewrite the first determining equation to become: 

-2h' + 2AcrV"^ In \u\ = -2{Xx + A2)(crV"^)^ In |n| - Iha'a'^ - 2m'cr"^'u~\ 

As h = f{x),a = cr{x), r = r{t, x), A, A2 £ M, the above relation is equivalent to the following ones: 

h' = a'a-^h, r = 0, Xa'a'^ = -{Xx + X2){a'a-^y . 

If a is an arbitrary function, then X = X2 = r = 0, h = Ca, C G M and we get {dt,cr{x)udu) as 
the maximal symmetry algebra. Hence, the extension of symmetry algebra is only possible when the 
function tp = ct'o""^ is a (non-vanishing identically) solution of equation 

(ax + + aV = 0, a,P e R, \a\ + |/?| / 0. 

If a 7^ 0, then at the expense of displacements by x we can get /? = 0, so that = mx~^, m 7^ 0. 
Integrating the remaining determining equations we get 

g = —2mx~^ \n\u\, f = mx~'^[mu In^ \u\ — (m — l)u In \u\ + nu], m 7^ 0, m,nG]R. 

The maximal invariance algebra of the obtained equation is the three-dimensional Lie algebra {dt, 
\x\"^udu,tdt + xdx) isomorphic to A3, 7. 

Next, if a = 0, then /? 7^ and ip = m, m 7^ 0. If this is the case, we have 

1 1 
q = In lul, f = — uln^ \u\ tiln Inl + nn, n ^ R. 

The maximal invariance algebra of the above equation reads as 

(9f,9a:,exp {^\^ '"'^«) 

and is isomorphic to ^43.2. 

Similarly we prove that the list of inequivalent equations of the form 1)2. 8|) admitting three- 
dimensional symmetry algebras is exhausted by the equations given below. Note that the presented 
algebras are maximal. This means, in particular, that maximal symmetry algebra of equation ()2.8() 
is, at most, three-dimensional. 

yl3.2-invariant equations 

1) Utt = Uxx +Ux\n\u\ + \u \t? \u\ — \uhi\u\+ nu (n G M) : {dt, S^:, exp {—\x) udu); 

2 

2) utt = Uxx + 'm.\hi \u\ — t]ux + ^u[(ln \u\ — t)(ln |n| — t — 1)] + nu{m > 0, 

n ^ R) : {dx,dt + udu,exp (— imx) udu)- 

A3,4-invariant equations 

1) Utt = Uxx + a;~"^[21n \u\ + mx^^t + n\Ux + x^^uln |n| 

+{mx~^t + n — 2)x~^n In + jui'^x'^t'^u + ^m(n — 3)x~^tu + px~'^u 
( m 7^ 0, n,p G M) : {tdt + xdx,x~^udu, dt — ^x~^ In \x\udu)- 



^3.5"iiivariant equations 

1) utt = Uxx + l^ip^ix + n|u|-^+^™' (m 7^ 0, n G M) : {dt, dx,tdt + xdx — m^^udu); 

2) Utt = Uxx + e'^Ux + ne^" (n G M) : ((9t, S^;, t^t + xdx - <9„); 

3) = Uxx — x~"^[21n |u| — mx~^t — n]ux + x~'^uln'^ \u\ 
-x~^{mx^^t + n)uln\u\ + ux~'^ [^x~'^t^ + -f (n - +p] 
{m,n,p £ M) : (t^t + xdx,xudu,dt + ^x'^udu). 

^3,7-invariant equations 

1) Utt = Uxx — 2mx^^Ux ln\u\ + mx^'^[mu\v? \u\ — [m — l)nln \u\ + nu] 

(m/0,l;neM) : {dt,\x\'^udu,tdt + xdx); 

2) tift = Uxx — x^^[2k + In \u\ — mx~^t — n]ux + k'^x^'^ulv? \u\ 

—kx~'^[mtx^^ + k + n — l]tiln \u\ + ^'m{k — 2 + n)tx~^u 
+ \mH'^x^^u + px~'^u {\k\ / 0, 1; m / 0, n,p G E) : 
{tdt + xdx, \x\''udu,dt + 2(i+fc) ^~'^^^«)- 

This completes group classification of nonlinear equations H2.8|) . 

4 Group classification of equation (2.12) 

Omitting calculation details we present below the determining equations for symmetry operators 
admitted by equation H2.12() . 

Assertion 3 The maximal invariance group of PDE is generated by the infinitesimal operator 

Q = T{t)dt + ax)dx + Ht)u + r{t, x)]du, (4.1) 
where r, ^, h, r, /, g are smooth functions satisfying the conditions 

rtx + f[h -n- ^x] = grx + Tft + ifx + [hu + r]/„, (4.2) 
ht = ng + rgt +igx- 

Assertion 4 The equivalence group £ of \2.1'J\) is formed by the following transformations of the 
space V; 

(1) t = T{t), x = X{x), V = U{t)u + Y{t,x), t'X'U^O; (4.3) 

(2) t = T{x), x = X{t), v = ^{x)^{t,x)u + Y{t,x), t'X'^^0, 

<i>{t,x) = exp[- J g{t,x)dt], 5^/0. 

As the direct verification shows, given arbitrary functions g and /, it follows from (|4.2|) that 
r = /i = ^ = r = 0. So that in the generic case the maximal invariance group of H2.12() is the trivial 
group of identical transformations. 

We begin classification of H2.12() by constructing equations that admit one-dimensional symmetry 
algebras. 



Lemma 3 There exist transformations i4-^ reducing operator i4.1\ ) to one of the seven canonical 
forms given below 

Q = tdt + xdx] Q = dt] Q = dx + tudu] 

Q = d^ + eudu, e = 0,l; Q = tud^, (4.4) 
Q = udu, Q = r{t,x)du, r ^ 0. 

Proof. Transformations (|4.3|1 reduce operator Q (|4.H) to become 

Q = TT'd-t + iX'ds + [{tU' + Uh)u + TYt + + Ur]d^. (4.5) 
Provided a • ^ 7^ 0, we can choose as T, X, [/, Y non vanishing identically solutions of equations 

tT' = T, iX' = X, TU' + hU = 0, TYt + ^Y^ + Ur = 

thus getting operator Q (|4.5|) in the form Q = td^ + xdx- If r 7^ 0, a ^ = 0, then taking as T, U, Y 
solutions of equations 

tT' = 1, tU' + hU = {U ^ 0), TYt + Ur = 

reduces operator 1)4. 1() to become Q = df. li t = 0, ^ / 0, then under /i' 7^ we get operator 
Q = dx + ivdv Next, if h' = 0, we arrive at the operator Q = dx + evdv, where either e = or e = 1. 

Finally, the case r = ^ = 0, gives rise to the following operators Q = tvdy, Q = vdy, Q = r(t, x)dv. 
After rewriting these in the initial variables we get the operators listed in the statement of lemma. 
The lemma is proved. 

Theorem 6 There exist, at most, three inequivalent nonlinear equations \2.1}3\) that admit one- dimen- 
sional invariance algebras. The form of functions f, g and the corresponding symmetry algebras are 
given below. 

{tdt + xdx) ■■ g = t~^g{uj), f = t~'^f{u,uj), w = tx"\ guj^O, fuu^O] 
(dt) ■■ g = g{x), f = f{x,u), g' 7^ 0, /„„ 7^ 0; 
{dx + tudu) :g = x + g{t), f = e^''f{t,Lo), uj = e'^^'u, f^^ ^ 0. 

Proof. If equation (|2.12|) admits one-parameter transformation group, then the latter is generated 
by infinitesimal operator H4.1() . According to Lemma 01 there exist equivalence transformations (|4.3() 
reducing this operator to one of the seven canonical operators 1)4. 4p . With this fact in hand, we turn 
to solving the determining equations (|4.2|1 for each of those operators. The first three operators yield 
invariant equations and corresponding symmetry algebras given in the statement of theorem. The 
next two operators give rise to inconsistent equations. 

Finally, the remaining operators yield that the functions / and g are linear in u, which means that 
the corresponding invariant equations are linear. 

It is straightforward to verify that for the case of arbitrary functions /, g, the corresponding 
one-dimensional algebras are maximal invariance algebras. 

The theorem is proved. 

We proceed now to analyzing equations (|2.12j) admitting two-dimensional symmetry algebras. 

Theorem 7 There exist, at most, three inequivalent nonlinear equations \2.1i^) that admit two-dimen- 
sional symmetry algebras, all of them being A2.2-invariant equations. The forms of functions f and g 



A\ = 
Al = 
Al = 



and the corresponding realizations of the Lie algebra A2,2 are given below 

^2.2 = i'^^t + xdx,t^dt + x^dx + mutdu) (m G M) : 
g=[mt + {k- m)x\t~^{t - A; / 0, 

/ = |t-xr-2|x|-"^/(a;), 
a; = n|t-x|-™|xr,/^^/0; 
= {tdt + xdx,t^dt + mtudu) (m G M) : 

g = t-'^[kx + mt], k^o, / = itr-vr'"/M, 

uj = \t\-"'\xru, L^^O; 
^2.2 = + xdx,x^dx + : 

g = (txy^{mx - t) (m G M), / = x"^ exp(-tx"^)/(w), 
a; = uexp{tx'^), f^^ / 0. 

To prove the theorem we need to extend reaUzations A\ (i = 1,2,3) to reahzations of the algebras 
^2.1,^2.2 by operators (|4.H) . We skip calculation details. 

Note that if the functions / are arbitrary, then the invariance algebras given in the statement of 
Theorem 13 are maximal. 

Now we can complete the group classification presented in Theorem |7| with the use of Lie-Ovsyan- 
nikov classification routine. 

We consider in some detail the case of A2 2-invariant equations (the remaining cases are treated in 
a similar way). The second determining equation from (|4.2j) reads now as 

{t - xfht = t''^Tt[m{t - x)^ + kx{t - x)] + T[-t~^m{t - xf - 2kt~^x + kf'^x'^] + k^,. (4.6) 

Differentiating right- and left-hand sides of (|4.(i|) twice by x yields 

ht = im-k){t-^Tt-t-\) + kf. 

Hence we get = and 

C = Xix'^ + X2X + X3, Ai,A2,A3GM, 

h = {m - k)t~^T + Xikt + X4,, A4 G M. 

With account of the above results we obtain from (|4.6|) that r = Ait^-|-A2t-|-A3. So it follows from (|4.6() 
that the coefficients of infinitesimal operator (|4.1j) . which generates symmetry group of A2 2-invariant 
equation, read necessarily as 

T = Xit^ + X2t + X3, 

= Xix"^ + X2X + X3, 
h = mXit + {m — k)X3t^^ + {m — k)X2 + X4, Ai, A2, A3, A4 G M. 

Consequently, the first determining equation from ()4.2() takes the form 

{x~^{t - x)"^[((m - k)X2 + Xi){tx - x^) - (m - k)X3t-'^x'^ 
-kX^x + mX3t]uj + r\t - x|-"|xr}/^ 

-x-^{t - x)~H((m - k)X2 + X4){tx - x^) (4.7) 

— (m — A;)A3t~^x^ — A:A3X -|- mX^t]/ 

= \t- xr™+2|^r[rt^ - t-\m + kx{t - x)-V^]. 



It follows from (|4.7I) that if the functions/ are arbitrary, then the maximal invariance algebra of the 
equation under study coincide with the realization A2 2- What is more, an extension of the invariance 
algebra is only possible when the function / obey the following equation: 

{auj + b)f^ - af = c, (4.8) 

where a,b,c G R, \a\ + |6| / 0. On the other hand, it follows from (|4.8j) 

(aw + b)f^uj = 0, 

whence fi_juj = 0. We arrive at the contradiction, which proves that there are no extension of the 
realization A2 2 in question to the higher dimensional invariance algebra of the equation H2.12() . Analy- 
zing A2 2" and A2 3-invariant equations we arrive at the same conclusion. 

Consequently, there are no nonlinear equations of the form (|2.12|) whose maximal invariance alge- 
bras are solvable Lie algebras of the dimension higher than two. Next, as the algebra s/(2,R) contains 
two-dimensional subalgebra isomorphic to A2.2, there are no nonlinear equations (|2.12|) . whose invari- 
ance algebras are either isomorphic to s/(2,M) or contain it as a subalgebra. Finally, we verified that 
there are no realizations of the algebra so(3) by operators (|4.1|) . 

Summing up the above reasonings we formulate the following assertion. 

Theorem 8 A nonlinear equation \2.1i^) having non-trivial symmetry properties is equivalent to one 
of the equations listed in Theorems and ^ 

This completes group classification of the class of nonlinear PDEs (|2.12|) . 

5 Group classification of equation (2.13) 

As earlier, we present the results of the first step of our group classification algorithm skipping deriva- 
tion details. 

Assertion 5 Invariance group of equation i2.1^) is generated by infinitesimal operator 

Q = T{t)dt + ^{x)d^ + {ku + r{t, x))du, (5.1) 

where k is a constant r, ^, r, / are functions satisfying the relation 

rtx + [k-T'- i']f = Tft + Ux + [ku + r]fu. (5.2) 

Assertion 6 Equivalence group E of the class of equations h2.1'J\) is formed by the following transfor- 
mations: 

(1) t = T{t), x = X{x), V = mu-hY{t,x), 

(2) t = T{x), x = X{t), V = mu + Y{t,x), T'X'm^O. (5.3) 

Note that given arbitrary /, it follows from (|5.2I) that T = ^ = A; = r = 0, i.e., the group admitted 
is trivial. To obtain equations with nontrivial symmetry we need to specify properly the function /. 
To this end we perform classification of equations under study admitting one-dimensional invariance 
algebras. 

Lemma 4 There exist transformations from the group E that reduce j|) to one of the four 

canonical forms: 

Q = dt + d^ + eudu (e = 0,1): 

Q = dt + eudu (e = 0, 1); 

Q = udu, Q = g{t,x)du {g j^O). 



Proof. Utilizing transformations (1) from (|5.3j) we reduce the operator Q to one of the foHowing forms: 



Q = di + ds + evd^ (e = 0, 1) : 

Q = dt + evd^ (e = 0,l); 

Q = ds + evd^ (e = 0, 1); 

Q = vdv, Q = g{t,x)dv {g 0). 

Next, we note that the change of variables t = x, x = t, v = v, which is of the form (2) from (|5.3jl . 
transforms the second operator into the third one. 

Finally, rewriting the obtained operators in the initial variables completes the proof. 

Theorem 9 There exist exactly two nonlinear equations of the form \2.1°J\) admitting one- dimensional 
invariance algebras. The corresponding expressions for function f and invariance algebras are given 
below. 

A\ = {dt + d^ + eudu) {e = 0,l): f = e''f{9,u;), = t - x, co = e-^'u; fl^ ^ 0; 
Al = {dt + eudu) (e = 0, 1) : / = e''f{x,uj), u = e'^^^x, /^^ / 0. 

To prove the theorem, it suffices to select those operators from the list of Lemma 0] that can be 
invariance algebra of nonlinear equation of the form 1)2. To this end we need to solve equation 
(|5.2|) for each of the operators in question. 

The first two operators yield A\- and Af-invariant equations. The last two operators gives rise to 
linear invariant equations (|2.13|) . which are not taken into account. 

What is more, if the function / is arbitrary, then the algebras A}^ and are maximal invariance 
algebras of the corresponding equations. 

Next, we classify nonlinear equations admitting symmetry algebras of the dimension higher than 
one. We begin by considering equations whose invariance algebras contain semi-simple subalgebras. It 
turns out, that the class of operators (|5.1|) contain no realizations of the algebra so(3). Furthermore, 
it contains four inequivalent realizations of the algebra sZ(2,R) given below. 

(1) {du\e^'du-\e-^'dt)- 

(2) {^dt,\e^\dt + du),-\e-^\dt-du)y, 

(3) {dt, \e^\dt + -^e-2*(ai - 

(4) {dt + a„ y^'dt + \e^^d.,, -\e-^'dt - ^e'^^^, + efe'^^ - e-2*]9„), e = 0, 1. 

Before analyzing s/(2, M)-invariant equations let us briefly review the group properties of the Li- 
ouville equation 

Ute = Ae", A/0. (5.4) 

It is a common knowledge that the maximal invariance group of this equation is the infinite-parameter 
group generated by the following infinitesimal operator |4Uj : 

Q = h{t)dt + g{x)d^-{h' + g')du, 

where h and g are arbitrary smooth functions. Note that due to this fact the Liouville equation can 
be linearized by a (non-local) change of variables (see, e.g., [T H 1 ^ Il2j). 



After a simple algebra we obtain that realizations (1), (3), (4) with e = 1 cannot be invariance 
algebras of nonlinear equation of the form H2.13() . Realization (2) is the invariance algebra of equation 



which reduces to equation (|5.4I) via the change of variables 

t = t, x = x, M = — -(f — In I/I), v = v{t,x). 

Finally making use of the change of variables 

t = e"^*, X = e"^"", V = u 

we rewrite (4) under e = to become 

{dt + d^,tdt + xd^,t^dt + x^d^). 

The corresponding invariant equation reads as 

ut^ = {t-x)-''j{u), fuu^O. (5.5) 

If the function / is arbitrary, then the above presented realization is the maximal invariance algebra of 
the equation under study. Using Lie-Ovsyannikov algorithm we establish that extension of symmetry 
is only possible when / = Ae" + 2. However, the corresponding equation is reduced to the Liouville 
equation by the change of variables 

t = t, X = X, u = v{t,x) + 2ln\t — x\. 

Thus the only inequivalent nonlinear equations (|2.13|) whose invariance algebras contain semi- 
simple subalgebras are given in 1)5. 4() and 1)5. 5() . where / is an arbitrary smooth function of u. 

To complete group classification of equation (|2.13|) we need to describe equations whose invariance 
algebras are solvable Lie algebras of the dimension higher than one. We begin with those realizations 
of two-dimensional Lie algebras A2.1, A2.2, which can be admitted by nonlinear equations (|2.13)) . 

It turns out that the class of operators (|5.1|) contains within the equivalence relation only one 
realization of the algebra A2.1 which meets the invariance requirements, namely, 

{dt + eiudu, + e2udu) (ei = 0, 1; €2 = 0, 1). 

The corresponding invariant equation reads as 

utx = ey^v{^i't + ^2x) f {oj) , u; = nexp(— eit — e2x). (5-6) 

Analysis of equation 1)5. 6|) with arbitrary /(w) shows that under ei -|- €2 7^ the above realization is 
its maximal invariance algebra. Provided ei = 62 = 0, the equation takes the form 

utx = f{u) (5.7) 

and its maximal invariance algebra is the three-dimensional Lie algebra of the operators 

{dt,dx,tdt - xdx), 

which is isomorphic to ^3.6- 

It is a common knowledge (see, e.g., |17| I1S( I19j ) that (|5.7|) admits higher symmetry if it is 
equivalent either to the Liouville equation ()5.4() . or to the nonlinear d'Alembert equation 

^Xte = A|nr+\ A/0, n/0,-1. (5.8) 



The maximal invariance algebra of H5.8jl is the four-dimensional Lie algebra of the operators 



{tdt - -udu,xdx - -udu, dt,dx). 
n n 

It is isomorphic to the Lie algebra A2.2 © ^2.2- 

Extension of symmetry algebra of equation H5.6|) with ei = 1, £2 = 0, is only possible when: 



Utx 
Utx 



-mt\ \m+l 
I "I ) 



Ae exp 

The maximal invariance algebra of H5.9() is the four-dimensional Lie algebra of operators 

{dt + udu, e^^dt, dx,xdx - —udu), 

m 

which is isomorphic to ^42. 2 © ^2.2- Note that the change of variables 



(5.9) 
(5.10) 



-mt 



X = X, U = v{t, x) 



reduces the above equation to the form l|5.8|) . 

The maximal invariance algebra of ()5.1U|) is spanned by the operators 

{dt + udu, dx,xdx - e^du), 

and is isomorphic to Ai © ^2.2- 

Analysis of 742.2-invariant equations yields the following results. The class of operators ()5.ip 
contains six inequivalent realizations of the algebra ^2.2 which meet the invariance requirements 



(1) 
(2) 
(3) 
(4) 
(5) 
(6) 



-tdt + xdu,dt); 

-tdt - xdx,dt + dx); 

-tdt - xdx + udu, dt + dx); 

-tdt + du,dt); 

-tdt - xdx - udu,dt); 

-tdt -xdx,dt). 



(5.11) 



(5.12) 



Equation invariant under realization (1) reads as 

Utx = exp(x"^u). 

Its maximal symmetry algebra is the three-dimensional Lie algebra of operators 

{-tdt + xdu, dt, xdx + udu) 
isomorphic to A2.2 © Ai. Note that the change of variables 

i=x, x = e*, u = v{t,x) 

reduces 1)5. 12|) to the form (|5.1U|) . 

Equation invariant under the second realization of ^2.2 is of the form (|5.5|) . It has already been 
studied while describing sZ(2, M)-invariant equations. 

Realizations (3) and (4) give no new invariant equations as well. 

New invariant equation are obtained with the use of the fifth realization from (|5.1()|) . It has the 
form 

Utx = x''^f{uj), OJ = x~^u. 



If the function / is arbitrary, then the reahzation in question is maximal invariance algebra of the 
above equation. Further extension of symmetry properties is only possible if f{uj) = X\uj\"^~^^, which 
gives the following invariant equation: 

utx = X\x\-'^-^\ur+\ A^O, 0,-1, -2. 

Its maximal symmetry algebra is the three-dimensional Lie algebra having the basis 

Tn -\- 1 

{dt,tdt + xdx + udu, xdx H udu) ■ 

m 

This algebra is isomorphic to A2.2 © Ai. 

We sum up the above results in the following assertion. 

Theorem 10 The Liouville equation Utx = Ae", A 7^ 0, has the highest symmetry among equations 
i2.1cl\) . Its maximal invariance algebra is infinite- dimensional and is spanned by the following infinite 
set of basis operators: 

Q = h{t)dt + g{x)dx - {h'{t) + g'{x))du, 

where h and g are arbitrary smooth functions. Next, there exist exactly nine inequivalent equations of 
the form \2.1'-^) . whose maximal invariance algebras have dimension higher that one. Those equations 
and their invariance algebras are given in Table 1. 

Table I. Invariant equations 116. y\) 



Number 


Function / 


Symmetry operators 


Invariance algebra 
type 


1 


e*/M, 
u; = ue~\fu,uj / 


dt + udu,dx 


A2.1 


2 


UJ = ue-*-'', f^uj + 


dt + udu, 
dx + udu 


A2.1 


3 


iO=(t- X)U, fujcu / 


-tdt - xdx + udu, 
dt + dx 


A2.2 


4 


UJ = x~^u, f^^ / 


-tdt - xdx - udu, 
dt 


A2.2 


5 


(t-x)-^/(n), 

/«« /o 


dt + dx, 
tdt + xdx, 
t^dt + x^dx 


sl{2,R) 


6 


exp(a;^^n) 


-tdt + xdu, 
dt,xdx + udu 


A2.2 © Ai 


7 


X\x\-"'-'^\u\"'+', 
A / 0,m / 0,-,l - 2 


dt,tdt - ^udu, 
xdx + '^^^ udu 


A2.2 © Ai 


8 


f{u)Juu^O 


dt,dx, -tdt - xdx 


^3.6 


9 


A|n|"+i,A / 0,n / 0,-1 


tdt - ^udu 
xdx - ^udu 
dt,dx 


A2.2 © A2.2 



6 Group classification of equation (2.7) 

The first step of the algorithm of group classification of 1)2. 7() 

Utt = Uxx + F{t , X ,U,Ux) , Fu^Uj: 7^ 



has been partially performed in the second chapter. It follows from Theorem ^ that the invariance 
group of equation (|2.7() is generated by infinitesimal operator H2.4|) . What is more, the real constants 
A, Al, A2 and real-valued functions h = h{x),r = r{t, x), F = F{t, x, u, Ux) obey the relation (|2.5j) . The 
equivalence group of the class of equations 1)2 .Tf) is formed by transformations p.lf) . 

The above enumerated facts enable using results of group classification of equation (|2.8|) in order 
to classify equation (|2.7() . In particular, using Lemmas ^ and [21 it is straightforward to verify that the 
following assertions hold true. 

Theorem 11 There are, at most, seven inequivalent classes of nonlinear equations \2. 1\) invariant 
under the one- dimensional Lie algebras. 

Below we give the full list of the invariant equations and the corresponding invariance algebras. 

{tdt + xdx) ■ F = t~^G{£,,u,uj), ^ = tx^^ , oj = xux] 
{dt + kdx) (fc > 0) : F = G{r], u, Ux), r] = x - kt; 
(dx) : F = G{t,u,Ux); 
(dt) : F = G{x,u,Ux); 
{dt + fix)udu) (/ / 0) : 

F = -tf'u + t'^if'fu - 2tf'ux + e*^ G{x, v, w), 
V = e~*'^u, uj = u^^Ux — f' f^^ In 

(/ / 0) : F=-f-'f"uln\u\ 
-2f^^f'ux ln\u\ + f-^{f'fuln^ \u\ + uG{t,x,uj), 
oj = u^^Ux — f'f~^ In ; 

{f{t,x)du) (//O) : F = f-\fu-fxx)u + G{t,x,uj), 

U; = Ux- f'^fxU. 

Note that if the functions F and G are arbitrary, then the given algebras are maximal (in Lie's sense) 
symmetry algebras of the respective equations. 

Theorem 12 An equation of the form \2. 7| ) cannot admit Lie algebra which has a subalgebra having 
nontrivial Levi factor. 

With account of the above facts we conclude that nonlinear equations 1)2. 7|) admit a symmetry 
algebra of the dimension higher than one only if the latter is a solvable real Lie algebra. That is 
why, we turn to classifying equations (|2.7|) whose invariance algebras are two-dimensional solvable Lie 
algebras. 

As the calculations are similar to those performed in the third section, we present the final result 
only. Namely, we give the form of invariant equations and the corresponding realizations of the 
two-dimensional invariance algebras. 
I. 742.1-invariant equations 

{tdt + xdx, udu) : F = x^'^uG{(,,lo), 

^ = tX~^ , UJ = U~^XUx] 

{tdt + xdx, cj{Odu) (a / 0, e = tx-^) ■■ 
F = - - 2ia')u + G(e,u;)], 

UJ = S^a'u + axux] 



A\ 
Al 
Al 
Af 
Al 



Al 



Al 



^2.1 — 
^2.1 — 



All = {dt + kdx, udu) (fc > 0) : F = uG{v, oo), 

rj = X — kt, CO = u~^Ux', 
^2.1 = {dt + kdx, v{ri)du) {k > 0, t] = x - kt, / 0) : 

F = {k^ — l)(p"(p~^u + G{ri, uj), ui = (pUx — <p'u; 
Ail = {dt + kdx, dx + udu) (A; > 0) : 

F = G{u!, v), 1] = X — kt, uj = ue~^, v = u~^Ux', 
A-ii = {dt, dx) : F = G{u,Ux); 
Al.i = {dx, udu) ■ F = uG(t,u;), LO = u~^Ux; 
Ail = {dx, mdu) {'P^O): 

F = (p~^ip"u + G{t,Ux); 
Ail = {dt, du) : F = G{x,Uxy, 
Al'i = {dt, f{x)udu) (/ 7^ 0) : 

F = -u~^ul + uG{x,Lo)- 

LO = u~^Ux — f'f~^ In l^xl; 
Ah\ = {dt + fix)ud^, g{x)udu) (S = f-'f - g-^g' / 0) : 

F = —g~^g"uln \ u\ — 2g^^g'ux In \ u\ 

+g''^{g')^uln'^ \u\ - 2f6tUx + 2f6g' g'^^tulnlul 
+f6H\ + f{g-^g" - f-'ntu + uG{x,Lo), 

CO = u~^Ux — g'g~^ In \u\ — tf6; 
Ah\ = {dt + fix)udu, e'fdu) (/ / 0) : 

F = [f- tf" + t\ff]u - 2tf'ux + e'^Gix, u), 

LO = e^*f{ux - tf'u); 
Al\ = {f{x)udu, g{x)udu) {S = f'g - g'f + 0) : 

F = -u~^ul - S~^S'ux 

+S-'[f"g' - //>ln|^.| + uG{t,xy, 
A^A = {V{t)du, mdu) (<^V - V^' + 0) : 

F = ip-^ip"u + G{t, X, Ux), (/c'V - Vi^" = 0. 

II. A2.2-invariant equations 

^2.2 = {tdt + xdx, xudy) ■ F = x'''^uIil'^ \u\ 

—2x~^Ux In |n| + f^uGi^, lo), ^ = tx"^; 

LO = xu~^Ux — In 
AI2 = {tdt + xdx, MOdn) (ifiy^O, ^ = tx-^) : 

F = t-\l - e)v>-'m^' + ^^'> + t-^Gi^, uj), 

LO = xifUx + ^^'u; 



{di + kdx, exp{k~^ x)udu) (/c > 0) : 

F = k~'^u\n^ \u\ - 2k~^Uxln\u\ - k~'^u\n\u\ 

+uG{ri,uj), 7] = X — kt, UJ = vr^Ux — k~^ In \u\; 

{dt + kdx, e*ip{ri)du) {r] = x - kt, k > 0, / 0) : 

F = {{k^ - l)<p"<p-^ - 2k(p'(p-^ + l)u + G{r],Lo), 

I I dtp 
u = ipux - ^pu, ip = — ; 

ar) 

i-tdt-xdx, dt + kdx) {k>0): 

F = r]~'^G{u,uj), 7] = X — kt, u = UxTj; 

{—tdt — xdx + mudu, dt + kdx) (fc > 0, m / 0) : 

F = \r]\-^-'^G{v,uj), ri = x-kt, 

u; = u\vr, v = ux\vr+'; 

{dx, e^udu) -F = nln^ |u| — "uln |m| — 2na; In l^^l 

+uG{t, uj), UJ = u~^Ux — In \u\; 

{dx, eXi)5„) (ip^O): 

F = {(p'''^(p" — l)u + G{t, uj), uj = Ux — u; 

{-tdt - xdx, dx) : F = t~^G{u,tUx); 

{-tdt - xdx + kudu, dx), (kj^O): 

F = \t\-^-''Giv,uj), V = \t\^u, UJ = \t\''+^ux; 

{dt, e^du) : F = u + G{x,Ux); 

{—tdt — xdx, dt) : F = x^'^G{u,uj), uj = xux] 

{dt + fix)udu, e^'+f^'du) (/ / 0) : 

F = - {tf - t\f)^ - (1 + f)) u - 2tfux 

+e'fG{x,u:), UJ = e-'f (n, - f{t + f-^)u) ; 

{—tdt — xdx, dt + kx^^udu) (A; > 0); 

F = -2ktx^^u + kH'^x~\ + 2ktx~'^Ux 

+x~'^ exp{ktx~^)G{v,uj),v = exp{—kx~^t)u, 

UJ = xu~^Ux + In |n|; 

{k{tdt + xdx), Ixl'''' udu) {k ^0,1): 

F = -k~^{l - k)x'~^u In |u| - 2k~^x~^Ux In \u\ 

+k~'^x~'^u\T? \u\ + x~^uG{v,uj), 
V = tx~^, UJ = xu~^Ux — k^^ In \u\; 

{k{tdt + xdx), \t\'''\{i)du) (^/o, 

i = tx-^) : F = [k-^{k-^ - 1) + 2i{k-^ - C^)^"V' 
- V"]i-'« + i-'G(C, uj), 

UJ = XipUx + ^f'u. 



In the above formulas G stands for an arbitrary smooth function. As usual, prime denotes the 
derivative of a function of one variable. 



6.1 Group classification of equation 

Utt = Uxx — u^^u^ + A{x)u,j. + B{x)u In |m| + uD{t, x) 

Before analyzing equations H2.7() admitting algebras of the dimension higher than two we perform 
group classification of the equation 



Here A{x), B{x), D(t, x) are arbitrary smooth functions. Note that the above class of PDEs contains 
j42\-invariant equation. Importantly, class 1)6.11) contains a major part of equations of the form (|2.7() . 
whose maximal symmetry algebras have dimension three or four. This fact is used to simplify group 
classification of equations 1)2. 7|) . 

Lemma 5 If A, B and D are arbitrary, then the maximal invariance algebra of PDE h6.1]) is the 
two-dimensional Lie algebra equivalent to A^^ and \6.1\) reduces to A^^-invariant equation. Next, if 
the maximal symmetry algebra of an equation of the form \6.1\) is three-dimensional (we denote it as 
^3 ), then this equation is equivalent to one of the following ones: 



I. A^ ~ A3.1, A3 = {dt, f{x)udu,ip{x)udu), 

A = -a-^a', B = a'^p, D = 0,a = fip- ftp' ^ 0, 
p = ^'f"-^"f; 

II. A3 ~ yls.i, A3 = {f{x)udu, ip{x)udu, dt + ip{x)udu), 
A = -a'^a', B = a'^p, 
D = ta-^[a'^' -il)p-ail)"l 
a = f^-ip'f^O,p = r^'-^"r, 

///. D = x-^G{^), C = tx-^, G^O: 

1) A3 ~ A3.2, A3 = {tdt + xdx,udu, |x|^""n(9„), 

A = nx-^ (n / I), B = 0; 

2) A3~A3.3, A3 = {tdt + xd^,udu,uln\x\du), A = x-\ B = 0; 

3) A3 ~ A3.4, A3 = {tdt + xdx, ^/\x\udu, \x\udu), 

A = 0, i?=ix-2; 

4) A3 ~ A3.9, A3 = {tdt + xdx, y/\x\cos{\(j\n\x\)udu, 



sin(2/51n |x|)n(9t(), A = 0, i? = mx ^, 
m> \, (3 = V4m - 1; 

5) A3 ~ A3.7, A3 = {tdt + xdx, i^iY+^udu, {y^\y-Pudu), 

A = 0,B = mx~^, m < J, m / 0, (3 = — 4m; 

6) A3 ~ A3,8,A3 = {tdt x^x,cos{^/rnln\x\)u^u, 

sm{^/mln\x\)u^u), A = x~^ , B = mx~^ , m > 0; 

7) A3 ~ A3.6,^3 = {tdt + xdx, \x\^\udu, \x\-^\udu), 
A = x^^, B = mx^^, m < 0; 



Utt = Uxx — u ^u^ + A{x)ux + B{x)u In \u\ + uD{t, x). 



(6.1) 



8) As ~ ^3.4, As = {tdt + xd^, {y^\y-^udu, ( v1^)^-" 
X In \x\udu), A = nx~^ (n / 0, 1), B = \{n - 

9) As ~ ^3.9,^3 = {tdt + xd,,{^\f-^cos{\l3\n\x\)udu, 
{^\f-'^sin{\f3\n.\x\)udu), A = nx'^ (717^ 0,1), 

B = mx"^ (m > i(n - 1)^), (3 = sj Am - (n - 1)^; 

jo; A3 ~ A3.7, ^3 = (t^t + xd,, [^\f-^--ud^, ( vH)^-"+^ 

xudu), A = nx~^ (n 7^ 0, 1), B = mx~'^ 
[m < l{n - 1)^, m / 0), /3 = v^(n - 1)^ - 4m. 
/F. D = G{t), 

1) As ~ A3.3, ^3 = {dx,udu,xudu), 
A = B = 0; 

2) As = As.2, As = (a^,u5„,e^u5„), 

A = -l, 5 = 0; 

3) As ~ A3.8, As = {dx,cos{x)udu,sm{x)udu), 

A = 0, B = l; 

4) As ~ Asm, As = {dx,e^udu,e-''udu), 
A = 0, B = -l; 

5) As ~ Asa, As = (9a;,exp [\x) udu,ex.p [^x) xudu), 
A = -l, B = \; 

6) As ~ As.7, As = (5a;,exp (i(l + udu,ew - f3)x) udu), 

A = -l,B = m {m<l), m ^ (3 = ^/l - Am\ 

7) As ~ A3.9, As = {dx,ex.p (^x) cos{^Px)udu,exp [^x) sin(i/3x)u(?„), 

A = -1,B = m {m > I), P = ^/4m - 1; 
F. D = G{r]), r] = x-kt, k > 0, 

1) As ~ A3.3, As = {dt + kdx, udu, xudu), 
A = B = 0; 

2) As = As.2, As = {dt + kdj:, udu, e^'udu), 
A = -l, 5 = 0; 

3) As ~ As.8, As = {dt + kdx, cos {x)udu, sin {x)udu), 
^ = 0, 5 = 1; 

A3 ~ A3.6, A3 = {dt + A;5j;, e^n5„, e~^udu), 

A = n, B = -l; 
5) As ~ A3.4, A3 = {dt + kdx,ex-p [^x) udu,exp [^x) xudu), 

A = -l, B = ^; 



6) ^3 ~ -43.7, M = {dt + kd^,ew (^(1 + P)x) U(9„,exp (i(l - udu), 
A = -1,B = m (m < I), m / 0,/3 = ^/T^^lm; 

7) A3 ~ ^3.9, A3 = {dt + A:5^.,exp (ij;) cos(i/3x)n5„, 

exp (ix) sin(^/3x)ii5u), A = —1, B = m {m > \) [3 = \J Am — 1 . 

Proof. Inserting the expression 

F = —u^^u^ + A{x)ux + B{x)u In |u| + uD{t, x) 

into classifying equation ()2.5p we get the system of determining equations for the functions h{x), r{t, x) 
and constants A,Ai,A2: 

r = 0, (Ax + X2)A' + XA = 0, 

(6.2) 

(Ax + X2)B' + 2XB = 0, h" + Ah' + Bh = -{Xt + Ai) A - (Ax + As)!?^ - 2AZ?. 

First, consider the case of arbitrary functions A,B,D. The left-hand side of the fourth equation 
of (|6.2j) depends on x only. What is more, since D is arbitrary, relation Dt ^ holds. Hence it 
immediately follows that the constants A,Ai,A2 must be equal to zero. As a consequence, the fourth 
equation becomes linear ordinary differential equation for the function h{x) 

h" + Ah' + Bh = 0. (6.3) 

The general solution of the above equation reads as 

h = Ci/(x) + C72V9(x), Ci, C2 G E, 

/(x) and <f{x) being the fundamental system of solutions of the equation 

y" + Ay' + By = 0, y = y{x). (6.4) 

Inserting this expression into 1)6. 3() yields 

A = -a-'a', B = a-\^'f"-f'^"l 

where a = iff — (p' f 7^ 0, which proves the first part of lemma. 

Suppose now that D = Q. Then, if at least one of the functions A ox B is arbitrary, then A = A2 = 
and the function /i is a solution of H6.3|) . This completes the proof of the case I of the second part of 
the lemma statement. 

Provided functions A and B are not arbitrary, it follows from the second and third equations of 
(|6.2|) that one of the following relations 

1) A = B = {); 

2) j4 = n, i? = m, m, n S M, |n| + \m\ 7^ 0; (6.5) 

3) A = nx^^, B = mx~'^, m,n S M, \n\ + \m\ ^ 

holds. With these conditions the maximal invariance algebra of (|6.1j) has the dimension higher than 
three. Consequently, without loss of generality we can suggest that D ^ 0. Integrating the equation 

{Xt + Ai)A + (Ax + X2)Dx + 2XD = H{x), 



under D ^ yields the following (inequivalent) expressions for the function D(t,x) 



D = x-'^G{C) +x-'^ j xH{x)dx, i = tx-^] 

D = G{'n) + k-^ fH{x)dx,'n = x-kt, k>0; 

D = G{t) + J H{x)dx, ^^'^^ 

D = tH{x) + H{x). 

The change of variables 

t = t, x = x, u = e{x)v(t,x), e^O, (6.7) 
where is a solution of equation 

0-^0" _ 0-2(0')2 + j^g-i0' ^Bln\e\+ A(x) = 0, 

preserves the form of equation (|6.1() . We can use this fact to simplify the form of the function D. As 
a result, we get 

D = e = tx-^; 

D = Girt), r] = x-kt, k > 0; 

(6.8) 

D = G{t), 
D = tH{x). 

If the function D is given by the one of the first three expressions, then H{x) = and h satisfies 

Given the condition D = tH{x), we have 

h" + Ah' + Bh = -XiH, {Xx + X2)H' + 3XH = 0. 

So that the maximal invariance algebra of the corresponding equation (|6.1() is three-dimensional iff 
A = A2 = 0, which yields the case II of the second part of the lemma statement. 

Turn now to the case when D = x~^G(^), ^ = tx~^. Then the function G ^ obey the equation 

{X2^ - Xi)G' + 2X2G = 0. (6.9) 

If G is an arbitrary function, then Ai = A2 = 0. In addition, we have A 7^ (otherwise the maximal 
invariance algebra is two-dimensional). Hence we get 

xA' + a = 0, xB' + 2B = 0. 

Consequently, functions A and B are given by either first or third formula from (|6.5j) . Analyzing these 
expressions yields ten cases of the case III of the second part of the lemma statement. 
If the function G is not arbitrary, then integrating (|6.9j) we get 

G = p, peR, py^O; 

G = p{i-q)-\ p^O, (7>0. 

Given the condition G = p, the parameter A2 vanishes. Hence in view of the requirement for the 
maximal algebra to be three-dimensional, it follows that A vanishes as well. This yields the case when 
A and B in (|6.H) are arbitrary functions (the case I of the second part of the lemma statement). If 
G = — q)^"^, P 7^ 0, then Ai = X2q- Hence we conclude that the maximal invariance algebra of 
the corresponding equation 1)6. 1(1 is three-dimensional iff the functions A, B are given by formulas 3) 
from (|6.5)) (which implies that Ai = A2 = 0) and we get the case III of the second part of the lemma 



statement. Next, if A, B are given by formulas 2 from (|6.5|) (which imphes that A = 0, D = p{t—qx)^^) 
and we arrive at the case IV {q = 0) or the case V (g > 0) of the second part of the lemma statement. 
Turn now to the case D = G{r]), r] = x — kt, k > 0. If these relations hold, then 

X{riG' + 2G) + (Aa - kXi)G' = 0. 

Hence it follows that if G is an arbitrary function of rj, then A = 0, A2 = kXi. That is why, the 
maximal invariance algebra of ()6.1() is three-dimensional iff either A = B = or A,B are given by 
formulas 2 from 1)6. 5|1 . So that we have obtained all equations listed in the case V of the second part 
of the lemma statement. 

The cases when either G = p {p ^ 0) or G = prj^'^ {p ^ 0) yield no new invariant equations ()6.5() . 
Consider now the last possible case D = G{t). If this is the case, then the equation 

{Xt + Xi)G' + 2XG = 

holds. Hence if follows that if G is an arbitrary function, then A = Ai = 0. So that the maximal 
invariance algebra of equation 1)6. 1|) is three-dimensional iff A, B are given by the formula 2 from (|6.5() 
and we have obtained all the invariant equations from the case IV of the second part of assertion of 
the lemma. If either of relations G = p {p ^ 0) ov G = pt~'^ {p 7^ 0) hold, them no new invariant 
equations having three-dimensional maximal invariance algebras can be obtained. 

To complete the proof of the lemma, we need to establish non equivalence of the obtained invariant 
equations. To this end it suffices to prove that there are no transformations from the group £, reducing 
their invariance algebras one into another. 

As we already mentioned in Section 3 there exist nine non-isomorphic three-dimensional solvable 
Lie algebras A^i = (61,62,63) (i = 1,2,... ,9). We analyze in some detail the case of the algebra 
j43.3. The list of invariant equations and algebras contains three algebras which are isomorphic to 
^3.3, namely, 

Li = {tdt + xdx,udu,uln\x\du); 

L2 = {dx,udu,xudu); 

L2 = {dt + kdx,udu,xudu) {k > 0). 

Denote the basis elements of the algebra L2 as 61,62,63. Suppose that there is a transformation ip 
from the group £ transforming L2 into L3. In other words we suppose that there exist constants 
ai,f3i,6i & R {i = 1,2,3) such that the relations 



V'(ei) 



aiii, 99(62) = ^ Aei, ip{e3) = ^ SiCi 



i=l 



i=l 



i=l 



and 



A 



Pi 

Si 



Q2 03 

hold. In the above formulas ei = c% + kdx, 62 = vdy, 63 = xvdv Equating the coefficients of the 
linearly independent operators d^, dx, dy yields that ai = f3i = 61 = 0. Hence we get the contradictory 
equation A = 0. This means that realizations L2 and L3 are non-isomorphic. Analogously, we prove 
that Li and L2 (as well as Li and L3) are non isomorphic. 

The remaining algebras are considered in a similar way. The Lemma is proved. 

In what follows we will use the results on classification of abstract four-dimensional solvable real 
Lie algebras = (61,62,63,64) j43| I44j . There are ten decomposable 



4yli = 3^1 0^1 = ^3.1 ® Ai, A2.2 2^1 = yl2.2 © ^2.1 = ^3.2 © ^1, 
2^2.2 = ^2.2 © ^2.2, A3,i e Ai (i = 3, 4, . . . , 9); 



and ten non-decomposable four-dimensional solvable real Lie algebras (note that we give below non- 
zero commutation relations only). 



■ [62,64] =ei, [63,64] =62; 

^4.2 : [61,64] =gei, [62,64] =62, 

[63,64] = 62 + 63, g / 0; 

: [61,64] =61, [63,64] =62; 

A4.4 : [61,64] =61, [62,64] = 61 62, 

[63,64] = 62 +63; 

A4.5 : [61,64] =61, [62,64] = ge2, 

]63,64]=P63, -l<p<q<l, p-Qy^O; 

Ai,6 ■■ [61,64] = Q'ei, [62,64] =pe2 - 63, 

[63,64] = 62 +P63, Qy^O, p>0; 
M.l ■ [62,63] =61, [61,64] =261, 

[62,64] =62, [63,64] = 62 + 63; 

^4.8 : [62,63] =61, [61,64] = (1 q')6i, 

[62,64] = 62, [63,64] = 0^63, \q\ < 1; 

^4.9 : [62,63] =61, [61,64] = 2g6i, 

[62, 64] = 0^62 - 63, [63, 64] =62 + qe3, q > 0; 

Aa.IO ■ [61,63] =61, [62,63] =62, 

[61,64] = -62, [62,64] = 61. 

Theorem 13 Equation utt = Uxx — u~^u^ has the widest symmetry group amongst equations of the 
form \6.1\) . Its maximal invariance algebra is the five- dimensional Lie algebra 

A\ = {dt, dx,tdt + xdx, xudu,udu). 

There are no equations of the form \6. 1]) which are inequivalent to the above equation and admit 
invariance algebra of the dimension higher than four. Inequivalent equations admitting four- 

dimensional algebras are listed below together with their symmetry algebras. 

I. D = 0, 

1) Aa ~ ^3.6 ©^1, M = {dt, dx, uch{Px)du, usinh{(3x)du) , 
A = 0,B = -P^, /3^0; 

2) A4 ~ ^3.8 ©^1, A4, = {dt,dx,ucos{/3x)du,usm{f3x)du), 
A = 0,B = f3^, /3/0; 

3) ^4- ^2.1 ©^2.2, A4 = {dt,dx,udu,e-'=udu),A = l,B = 0; 

4) ^4-^3.4©^!, A4 = {dt,dx,e~^udu,xe-''udu), A = 2,B = 1; 

5) A4 ~ j43.9 (B Ai, A4 = {dt,dx,ue~^ cos{Px)du,ue~^ sm{Px)du), 
A = 2,B = m,m > 1, j3 = \Jm - 1; 

6) Ai ~ A3.7 ©^1, A4 = {dt, dx,ue~'^dLi{l3x)du,ue'''^ s'mh{(3x)du) , 
A = 2, B = m, m > 1, m / 0, /3 = \/l — m; 



1) A4^A4^,2, A4 = {dt,tdt + xdr„s/\x\udu,us/\x\ln\x\du), 
A = 0, B = |x-2; 

8) ^4~^4.5, A4 = {dt,tdt + xdx,\x\^+'^udu,\x\^''^udu), A = 0, 



B = mx ^, m < i, m / 0, /3 = — m; 
9) A4 ~ A4,Q, A4 = {dt,tdt + xdx, 

^/\x\cos{/3ln \x\)udu, -^/|x[sin(/?ln \x\)udu), 

A = 0, B = mx~'^, m> \, = 



10) ^4^^4.3, A4 = {dt,tdt + xdx,u\n\x\du,udy), A = x ^, B = 0; 

11) ^4^^3.7 0^1, A4 = {dutdt + xd^,\x\^-''udu,udu), 
A = nx-\ B = 0,n / 0, 1; 

12) ^4-^4.5, A = (5t,t5t + x9^,|a;|5(i-")^x5„,|x|^(i-")uln|x|5„), 
A = B = \{n- ifx''^, n / 0, 1; 

13) ^4-^4.5, ^4 = (5t,t5t + x5a.,|x|5(l-"+/^)ua„,|x|5(l-"-/^)u5„), 

A = nx~^, B = mx"'^, m < j{n — 1)^, m 7^ 0, n 7^ 0, 
/3= ^{n - 1)2 - 4m; 

^4-^4.6, A4 = (5t,t5t + x9,,|x|5(i-")cos(/31n|x|)tia„, 

2(1""-) sin(/31n |x|)7i9ij), y4 = nx~^, B = mx~'^, 
m 7^ 0, n / 0, m > ^(n - 1)^, /3 = -^m - ^{n - 1)2; 

D = ktx-^, k>0, 

1) A4 ~ ^4,1, A4 = ((9t — \kx~^udu,tdt + xdx,xudu,udu), A = B = 0; 

2) A4 ~ A4.2, A4 = {dt — ^kx^^udu, tdt + xdx, -^/|x[m5u, -^/|x[ln \x\udu), 
A = 0, B = ix-2; 

3) A4 ^ A4.5, A4 = {dt - ^^x~'^udu,tdt + xdx, \x\^^^udu, \x\^~'^udu), 
A = 0, B = mx"^, m 7^ 0, -2, m < \, P = ^J\-m; 

4) ^4 ~ -44.2, A4 = {dt + lkx-^{l + 3ln\x\u)du,tdt + xdx, 
x'^udu,x-^udu), A = 0,B = -2x~'^; 

5) A4 ~ A4.6, A4 = {dt 

~ r^\-2^ ^''^duttdt + xdx, \/\x\uQ.os{l3\n.\x\)du, 



m+2 

\x\usm.{(3\n\x\)du) , A = Q,B = mx~'^, m>\,(3—y,,o ^, 



6) A4 ~ ^4.3, A4 = {dt - kx ^udu, tdt + xdx, udu,ulii \x\du), 
A = x-^,B = 0; 

7) A4 ~ A3.4 e Ai, A4 = {dt + kx-\l + In |x|)u5„, tdt + xdx, ud^, 
x-^udu), A = 2x-^, B = 0; 

8) A4 ~ A3.7 e Ai, A4 = {dt + :^x-^udu,tdt + xdx,udu, |x|i-"n5„), 
A = nx-'^,B = 0, 0,1,2; 

9) A4 = A4A, A4 = {dt — \kx^^ \r? \x\udu, tdt + xdx, 
x~^udu,x~^\n.\x\udu), A = ?>x~^, B = x~^; 



10) ~ A4.2, A4 = {dt - j^x-^ud^,tdt + xd^, |x|2(i~")n9,, 
1x1^(1-") In A = nx"\ B = l)22;-2, n/0,3; 

|x| 2(i~"~^)u9„), A = nx^^, B = mx^'^, 

n / 0, 2, m 7^ n - 2, m < i(n - l)^, /? = y^(n - 1)2 - 4m; 

iSj A4~^4.2, A4 = (t^t + + 3^x"Mn|x|tia„,x"^n(9u,|xp-"ua„), 

A = B = {n- 2)x-^ n / 0, 2, 3; 

^4 ~ ^4.6, ^4 = (i^t + xd^, dt - 2r^x~^udu, 
|x| 2(i~")ucos(/?ln |x|)c?„, |x| 2(i~"-)nsm(/31n |x|)(9„), 
A = nx^^, B = mx~'^, n 7^ 0, m 7^ 0, m > ;j(n — 1)^, 

///. D = kt, k>0, 

1) A4 ~ A4.1, A4 = - \kx'^udu, xudu,udu) , A = B = 

2) Ai ~ A4.3, A4 = (5a;, 9i - kxudu, e'^'udu^udu) , A = 1, S = 0; 

5^ ^4~A3.8©Ai, A4 = - A;/3~2n5„,'u,cos(/3a;)(9„,usin(/3a;)5„), 

^ = 0, 5 = /3 / 0; 
4) ^4~A3.6©^i, A4 = (Sa;,^* + A;/3~2u5„,'ach(/3a;)5«,nsmh(/3x)5„), 

^ = 0, 5 = -/32, /? / 0; 

5j ^4 ~ ^43.4 © Ai, A4 = {dx,dt — 4kudu,exp {—^x) udu,xexp {—^x^ udu), 
A = l, B = l; 

6) Ai ~yl3.7©^i, A4 = {dx,dt - km-^udu,eyi.Y> - I3)x) udu,ey.Y> + I3)x) udu) , 

A = l, B = m, m < I, m / 0, /3 = Vl - 4m; 

7j A4 ~ j43.9 © j4i, A4 = {dx,dt — km'^udujexp {—^x) cos{(3x)udu, 
exp (— ^x) sm.{f3x)udu) , 

-4 = 1, B = m, m > i, /? = y^m - i; 

/F. D = fct^^ A; 7^ 0, 

A4 ~ A4.8 {q = -1), A4 = {dx,tdt + xdx,xudu, udu), A = B = 0; 

V. D = m{x - kty"^, k> 0, m 7^ 0, 

^4 ~ ^4.8 (9 = -1), A4 = {dt + kdx,tdt + xdx,xudu,udu), A = B = 0. 

Proof. According to Lemma to get the list of inequivalent equations of the form H6.1|) we need 
to analyze the cases when either D = or D does not vanish identically and is obtained through 
equations ()6.8() . 

If D = 0, then function h satisfies equation (|6.3|) and functions A, B are given by one of the formulas 
H6.5p . It follows from 1)6.2^ that the highest possible dimension of an invariance algebra admitted by 
(j6.H) equals to five. An equation admitting this algebra is equivalent to the following one: 

utt = Uxx - u~^ul. (6.10) 

This proves the first part of the assertion of theorem. 



The remaining expressions for the functions A, B from (|6.5|1 yield the fourteen items of the case I 
of the second assertion of theorem. 

The fourth expression for the function D given in within the equivalence relation (|(i.7j) boils 
down to either D = ktx~^, k > or D = kt, k > 0. Here the function h satisfy one of the equations 

h" + Ah' + Bh = -XiH, 

where H = kx^^ or H = k. Analyzing of the corresponding expressions for functions A, B ()6.5|) 
results in the cases listed in the cases II and III. 

Next, the function D given by the third formula from 1)6. 8|) simplifies to D = kt~'^ {k ^ 0), whence 
we get the results listed in the case IV of the second assertion of theorem. Similarly, the second 
expression for the function D gives rise to formulas of the case V. 

The first expression for D from (|6.8|) gives no new invariant equations. 

What is left is to prove that the so obtained invariant equations are inequivalent. We omit the 
proof of this fact. 

Theorem is proved. 

6.2 Nonlinear equations (12. 7p invariant under three-dimensional Lie algebras 

The class of PDEs (|2.7|) does not contain an equation whose invariance algebra is isomorphic to a Lie 
algebra with a non-trivial Levi ideal fTheoremll2|). That is why, to complete the second step of our 
classification algorithm it suffices to consider three-dimensional solvable real Lie algebras. We begin 
by considering two decomposable three-dimensional solvable Lie algebras. 

Note that while classifying invariant equations 1)2. 7(1 we skip those belonging to the class 1)6. 1(1 . 
since the latter has already been analyzed. 

6.2.1 Invariance under decomposable Lie algebras 

As A3 I = 3Ai = A2.1 © Ai, ^3,2 = A2.2 © ^1, to construct all realizations of Aa.i it suffices to 
compute all possible extensions of the (already known) realizations of the algebras A2.1 = (ei, 62) and 
A2.2 = (ei, 62). To this end we need to supplement the latter by a basis operator 63 of the form (|2.4j) 
in order to satisfy the commutation relations 

[ei, 63] = [62,63] =0. (6.11) 

What is more, to simplify the form of 63 we may use those transformations from £ that do not alter 
the remaining basis operators of the corresponding two-dimensional Lie algebras. 

We skip the full calculation details and give several examples illustrating the main calculation 
steps needed to extend A2.1 to a realization of ^3.1. 

Consider the realization A^i- Upon checking commutation relations (|6.11|) . where 63 is of form 
((231), we get 

Ai = A2 = r{t, x) = 0, h = k = const. 

Consequently, 63 is the linear combination of ei, 62, namely, 63 = Aei -|- A;e2, which is impossible by 
the assumption that the algebra under study is three-dimensional. Hence we conclude that the above 
realization of A2 ^ cannot be extended to a realization of the algebra A31. 

Turn now to the realization A2 i ■ Checking commutation relations (|H.1H) , where 63 is of form (|2.4j) 
yields the following realization of ^3.1: 

{tdt + xda;,a{Odu,'j{Odu), ^ = tx~^, 
where 7V — ja' ^ 0. However, the corresponding invariant equation (|2.7j) is linear. 



Finally, consider the realization A2 i ■ Inserting its basis operators and the operator 63 of the form 
H2.4|) into and solving the obtained equations gives the following realization of vls.i : 

{dt,dx,udu). 

Inserting the obtained coefficients for 63 into the classifying equation H2.5|) we get invariant equation 

Utt = Uxx + uG{u;), LO = U~^Ux, 

where (to ensure non-linearity) we need to have Gu)w 7^ 0. 

Similar analysis of the realizations (i = 4, 5, . . . , 12, 14) yields three new invariant equations. 
For two of thus obtained ^3.1-invariant equations the corresponding three-dimensional algebras are 
maximal. The other two may admit four-dimensional invariance algebras provided arbitrary elements 
are properly specified. 

Handling in a similar way the extensions of A2.2 up to realizations of ^3.2 gives ten inequiva- 
lent nonlinear equations whose maximal invariance algebras are realizations of the three-dimensional 
algebra A3. 2 and four inequivalent equations (|2.7|) admitting four-dimensional symmetry algebras. 

We perform analysis of equations admitting four-dimensional algebras in the next sub-section. 
Here we present the complete list of nonlinear equations (|2.7|) whose maximal symmetry algebras are 
realizations of three-dimensional Lie algebras A^^^i and ^43.2. 
A3. 1 -invariant equations 

A\ ^ = {dt,dx,udu) ■■ 

F = uG{uj), uj = u~^Ux; 
Al^ = {dx,mdu,mdu) : 

a = Tp'ip — Tpif' ^0, cr' = : 

F = ip~^ip"u + G{t,Ux). 

743.2-invariant equations 

^3.2 = {dt,dx,e''udu) : 

F = —u~^Ux — nln \u\ + uG{uj), 

uj = u^^Ux — In \u\ : 
^3.2 = {-'t^dt - xdx, dt + kdx,udu) (A; > 0) : 

F = ur]~^G{oj), r] = X — kt, 

UJ = riu~^Ux; 
^3.2 = {-tdt-xdx + mudu,dt + kdx,\r]\~"'du) 

(f] = X — kt, k = m = or k > 0, m £ R) : 

F = m{k^ - l)(m + l)r]-'^u + \r]\'^-'^G{oj), 

UJ = \riY^{mu + rjUx); 
^3.2 = {dx-.e'udu, dt + mudu) (m > 0) : 

F = —u~^u^ — Ux + uG{uj), 

UJ = u~^Ux — In \u\ + mt; 



^3.2 



A 



3.2 



^3.2 



A 



3.2 



^3.2 



/ilO 
^3.2 



(-t(9t - xd^,d^,udu) : 

F = ut^'^G^uj), uj = tu~^Ux] 

{—tdt — xdx,dt + kx^^udu,udu) (/c > 0) : 

F = 2ktx''^Ux - 2ktx'^u + k'^fx^'^u + x~^uG{u;), 

uj = xu~^Ux + ktx~^; 

{—tdt — xdx,dt + kx^^udu, ex.p{ktx~^)du) {k > 0) : 

F = 2ktx'^Ux + {k^fx''^ - 2ktx'^ + k'^x'^)u + 

+x~'^ exp{ktx^^)G{LL!), u> = exp{—ktx~^){xux + ktx~^u); 

{:^idt + kdx),e'+'''du,e^du) (k > 0,rj = x - kt) : 
Ik 

F = (A;^ — l)u + G(r7, a;), a; = u^' — 

(5t + f{x)udu, e(i+^(^«*5,, /(x)e^(^)*a„) : 

F = -{tf - t\f'f - (1 + ff)u - 2tfux + e'fGix,uj), 

•-tf{ux-f'it + r^)u), f" + 2f + f = 0, //O; 

{k{tdt + xdx),\t\>''"\C\^du,\^\^du) {k^O;l) : 



UJ 



1 - k 



1-^2 



4A;2 



t-\ + t-^Gi^,uj) 



k — 1 



XUx + 



k - 1 

'2k 



-u 



, i = tx-\ 



6.2.2 Invariance under non-decomposable three-dimensional solvable Lie algebras 

There exist seven non-decomposable three-dimensional solvable Lie algebras over the field of real 
numbers. All those algebras contain a subalgebra which is the two-dimensional Abelian ideal. Conse- 
quently, we can use the results of classification of yl2.i -invariant equations in order to describe equations 
admitting non-decomposable three-dimensional solvable real Lie algebras. We remind that equations 
of the form (|6.1() has already been analyzed and therefore are not considered in the sequel. 

As an example, we perform extension of the realization A^^-^ to all possible realizations of non- 
decomposable three-dimensional solvable real Lie algebras. The remaining realizations are handled in 
a similar way. 

It is straightforward to verify that transformations 



t = 7t + 71, 



e7x + 72, 



V = p{x)u + 0{x), 



(6.12) 



where 7,71,72 G K, 7 7^ 0, e = ±1, p ^ are equivalence transformations for the realization 
^2^1 ~ {dt, f{x)udu) (/ 7^ 0). That is why, we may use the above transformation in order to simplify 
the form of operator 63. As a result we get three inequivalent expressions for 63 

1 ) 63 = tdt + xdx + r{t, x)du {n / or r = 0); 

2 ) e'i = dx + r{t, x)du (rt 7^ or r = 0); 

3 ) 63 = r{t, x)du (rt / or r = 1). 



Let ei = df, = f{x)udu and 63 = tdt + xd^ + r{t, x)du, then 

[ei, 63] = 5t + rtdu, [62, es] = -xf'udu - rfdu- 

Analyzing commutation relations for the algebras A^^i (f = 3, 4, . . . , 9) we obtain that the necessary 
conditions for A^^^ to admit extension to a realization of ^3.5 read as r = 0, xf = — /, of A^ q as 
r = 0, xf = f, and of A3.7 as r = 0, xf = — (0 < |g| < 1). So gives rise to the following 
realizations: 

^3.5 : ei = dt, 62 = x^^udu, 63 = tdt + xd^; 

: ei = dt, 62 = xudu, 63 = tdt + xd^; 
A3.7 : ei = dt, 62 = 63 = tdt + aJ^x (0 < \q\ < 1). 

If 63 = dx + r{t, x)du, then 

[61, 63] = [62, 63] = -fudu - rfdu- 

Analyzing commutation relations for A^ ^ (z = 3, 4, . . . , 9) we come to conclusion that the realization 
A^i cannot be extended to a realization of the above three-dimensional Lie algebras. 
The same conclusion holds true when 63 = r{t, x)du {rt 7^ or r = 1). 

Let 61 = f{x)udu, 62 = dt- If 63 = tdt + xdx + r{t,x)du (rt / or r = 0), then it follows from 
commutation relations 

[61, 63] = -{rf + xfu)du, [62, 63] =dt + ndu 
that the only possible extension of the realization is the realization of ^3.5: 

{x~^udu,dt,tdt + xdx). 

This realization coincide within notation with the already obtained one. 
Next, if 63 = dx + r{t, x)du (r^ / or r = 0), then 

[61, 63] = -{fu + rf)du, [62, 63] = rtdu- 

Analyzing commutation relations for 743.^ (i = 3, 4, . . . , 9) we come to conclusion that the realization 
^2^1 cannot be extended to a realization of the above three-dimensional Lie algebras. 
The same conclusion holds true for the case 63 = r{t,x)du (rt / or r = 1). 

Summing up the above considerations we see that the realization Ag'i can be extended to the 
following realizations of non-decomposable three-dimensional solvable real Lie algebras: 





~ ^3.5, 




= {dt, 


x~^udu,tdt + xdx); 




~ ^3.6, 




= {dt, 


xudu,tdt + xdx); 




~ A3.7, 




= {dt, 


\x\~'^udu,tdt + xdx) (0 < 



Solving the corresponding classifying equations yields the following invariant equations: 

: utt = Uxx ~ U'~^Ux — 2x~'^uln\u\ + x~'^uG{uj), u = xu~^Ux + iT^\u\; 
: utt = Uxx — u~^u'^ + x~'^uG{u>), 

oj = xu^^Ux — In \u\; 
: Utt = Uxx — u~^Ux — q{q + l)x^'^uln\u\ + ux~'^G{lo), 

LU = xu^^Ux + ql^lul (0 < Iq'I < 1). 

Note that the algebras L^,!?, are maximal (in Lie's sense) invariance algebras of the corresponding 
equations. 



While classifying nonlinear equations invariant non-decomposable three-dimensional solvable Lie 
algebras we discovered equations whose maximal invariance algebras are four-dimensional. For exam- 
ple, after extending the realization ^ up to a realization of the algebra A3 3 we got the following 
realization of the latter: 

{du,dt,d^ + tdu)- 

The corresponding invariant equation (|2.7|) reads as uu = Uxx + G{ux). However, the maximal in- 
variance algebra of the above equation is the four-dimensional Lie algebra {dt,tdu,du,dx), which is a 
realization of 743.3© Ai. Note also that we have obtained the above invariant equation when classifying 
743.1-invariant equations. 

By the above reason, we give below only those nonlinear invariant equations whose maximal 
symmetry algebras are three-dimensional non-decomposable solvable real Lie algebras. 
^3 3-invariant equations 

^3.3 = {udu, dt + kdx, mdt + k'^xudu) (A; > 0, m / 0) : 
F = —u~^u^x + uG{uj), uj = X — kt + mk^u'^Ux] 

^3.3 = {udu, dx,mdt + xudu) (m > 0) : 

F = —u~^u^ + uG{uj), UJ = t — mu~^Ux; 

^3.3 = {\t\^du,-\t\hn\t\du,tdt + xdx + ^udu) ■■ 

F = —-t~'^u + u^G{£,,uj), ^ = tx^^ , UJ = xu^; 

^1.3 = {du,-tdu,dt + kdx) (fe>0); 
F = G{r], Ux), rj = X — kt. 

A3.4-invariant equations 

AIa = {\vr-^du,dt + kdx,tdt + xdx + {mu + t\r]r-^)du) 
(rj = X — kt, /c > 0, m 7^ 1) : 

F = {k^ - l)(m - l)(m - 2)j]-'^u - 2k{m - 1)7/™^^ In |7?| 
+ |7?r-2G(w), = [j^ux - (m - l)^z]|r?| — ; 
^1.4 = {du, -tdu, dt + kdx + udu) (A: > 0) : 

F = e^G{r],u}), r] = x — kt, u) = e~^Ux; 

11 3 
^3.4 = {\t\2du,-\t\2\n\t\du,tdt + xdx + -udu) : 

F = -^t^'^u + u^^G{i, w), i = tx'^, UJ = x~^u\\ 
^3.4 = {kx~^udu, dt — kx^^ In \x\udu, tdt + xdx) (A; > 0) : 
F = —'iktx~^u — 2x^^nln \u\ — u^^u^ + x^'^uG{uj), 
UJ = xu~^Ux + In \u\ + ktx~^; 



^3.4 = {ex.]){ktx ^)du,dt + kx ^udu,tdt + xd^ + {u + texp{ktx (/c > 0) : 

F = k'^x-\{t^ + x^) + 2x-^{ktx-^ + l)ua; 
+2k exp(/ctx~"^)x~^ In |x| + x~^ ex.\>{ktx^'^)G(uj), 

CO = exp{—ktx~^){ux + ktx~'^u). 

Aa.s-invariant equations 

4.5 = (l^r"^5n, dt + kd^,tdt + xd^ + mudu) {k>0, m / 1) 

F = {k^- l)(m - l)(m - 2)?7-2u + {rj^-^Gito), 

u = \r]\~"^[r)Ux — (m — 77 = x — 
4.5 = {dt,dx,tdt + xdx) : 

F = ulG{u)- 
4.5 = (^t' ^a;, + + mudy) (m / 0) : 

F = \u\^-^G{u:), u = \uxr\u\'^-"'; 
4.5 = {9t, dx, tdt + xdx + du) : 

F = e-2"G((j), io = e'^Ux] 
4.5 = i^t, x~^udu, tdt + xdx) ■ 

F = -u~^ul - 2x"^uln \u\ + x~\G{uj), 

oj = xu'^Ux + In |n|; 
4.5 — i^t + kx~^udu, e-xp{ktx~^)du, tdt + xdx + udu) [k > Q) : 

F = kx'^u[kt'^ - 2tx + kx'^] + 2ktx~'^Ux + x~^ exp{ktx~'^)G{uo), 

u = exp{—ktx~^){ux + ktx~^u); 
4.5 = {ipit)du, ^(t)du, dx + udu) ((^V - V'V'' 7^ 0) : 

F = ip~^ip"u + UxG{t,cj), 

CO = e~^Ux, ip"'ip - (pip" = 0. 
Aa.Q-invariant equations 

4.6 = {dt + kdx, |r?r+'5„, tdt + xdx + mudy) {k >(), m ^ -1) : 

F = m{k'^ - l)(m + l)??-^^ + \rir-^G{co), 

^ — \v\^~"^Wx ~ 'r}~^{rn + = x — kt; 

4.6 ~ {^t + mx~^udu, xudu, tdt + xdx) > 0) '■ 

F = -u~^ul - 2mtx~^u + x~'^uG{co), 

CO = xu'^Ux — In Inl + 2mtx~^\ 



^3.6 — {dt + kx ^udu,exp{ktx ^)du,tdt + xd^ — udy) {k > 0) : 

F = x~'^[k'^x^ - 2ktx + k^t^]u + 2ktx~^Ux + x"^ exjp{ktx~^)G{u;), 
u = e'x.p{—ktx~^){x^Ux + ktu); 

4.6 = dt + kd^) {k>0): 

F = u + G{rj, Ux), rj = X — kt; 

4.6 = {\t\~^du,\t\Wu,tdt + xda; + ^udu) ■■ 

F=^t-\ + \t\-'^G{^,u;), i = tx-\ LO = x-'ul. 
^3 7-invariant equations 

4.7 = {dt + kd,,,\vr^''du,tdt + xdx + mudu) 

{k>0,m^q, < \q\ < 1) : 

F={k'^- l)(m -q)im-q- l)r]-^u + \r]r-^G{uj), 

^ — Ivl^'^b^x — {fn — q)rf^u], r] = x — kt; 

4.7 = {dt + kx'^udu,exp{ktx~^)du,tdt + xdx + qudu) 
{k>0, < \q\ < 1) : 

F = [k^x-^ + k^x'H^ - 2ktx'^]u + 2ktx~^Ux 
+ exp(A:tx-i)G(w), 

Lv = \x\^~'^ exp{—ktx~^){ux + ktx~'^u); 
4.7 = {\t\'^''d^,\t\'--^^dM+xdx + {l + ^q)ud^) iq^O,±l): 

F=^q{q-2)r^u+\t\-2^i-^^G{^,u), 

Al y = (exp Q(g - 5„,exp Q(l - 5„,at + + ^(1 + 
(g7^0,±l;A;>0) : 

F=\{q- Ifu + exp Q(l + q)t^ G{'n,u), 
ri = x- kt, Lo = exp (^-^(1 + 9)*^ "x; 

Al^ = {dt + kx-\du,\x\-''udu,tdt + xdx) {k>0,qj^O,±l) : 

F = -u~^ul - q{q + l)x"^uln \u\ + k{q - l){q + 2)tx~^u 
+ux~'^G{u!), UJ = xu~^Ux + gin \u\ + k{l — q)tx~^. 

^3 g-invariant equations 

4.8 = (cos tdu, - sin tdu, dt + kd^) (A; > 0) : 
F = —u + G{r], Ux),r) = x — kt; 



^3.8 = {\t\^^ cos{ln\t\)du,-\t\^^ sm{ln\t\)du,tdt + xd^ + ^udu) ■■ 
F = -^r^u + \t\-lGiC,u;), 

^ = tX , W = |t| 2 Ux- 

yl3 9 -invariant equations 

ylg 9 = (sin tdu, cos tdu, dt + kdx + qudu) {k > 0,q > 0) : 

F = -u + e^^G{r],uj), r] = x - kt, uj = e~'^^Ux] 
^3.9 = (|i|^sin(ln|t|)(9„,|t|3cos(ln|t|)5„,t9t + x5^ + (^ + g)u5„) 

(g/0):F = -V2n + |t|9-iG(e,a;), 
6.3 Complete group classification of equation (12.71) 

The aim of this subsection is finalizing group classification of H2.7() . The majority of invariant equations 
obtained in the previous subsection contain arbitrary functions of one variable. So that we can utilize 
the standard Lie-Ovsyannikov approach in order to complete their group classification. 

6.3.1 Equations depending on an arbitrary function of one variable. 

Note that equations belonging to the already investigated class of (|6.1|) are not considered. 
As our computations show, new results could be obtained for the equations 

Utt = Uxx + uG{uj), uj = u^^Ux, (6.13) 
Utt = Uxx + G{ux) (6.14) 

only. Below we give (without proof) the assertions describing their group properties. 

Assertion 7 Equation i6.1S\) admits wider symmetry group iff it is equivalent to the following equa- 
tion 

Utt = Uxx + m-u'^Ux (^1^/0,— 1). (6.15) 
The maximal invariance algebra of i6.15\) is the four- dimensional Lie algebra 

A4 ~ A3.5 © ^1, ^4 = {dt, dx,tdt + xdx,udu). 

Assertion 8 Equation ^6.14^ admits wider symmetry group iff it is equivalent to one of the following 
PDEs: 

Utt = Uxx + e^""; (6.16) 
Utt = Uxx + mlnluxl, m>0; (6.17) 
Utt = Uxx + \ux\\ fc/0,1. (6.18) 

The maximal invariance algebras of the above equations are five- dimensional solvable Lie algebras listed 
below. 

Al = {dt,dx,du,tdu,tdt + xdx + {u-x)du); 

Al = {dt,dx,du,tdu,tdt + xdx + {2u + -mt'^)du); 

k — 2 

A^ = {dt,dx,du,tdu,tdt + xdx + -^—judu). 



Analyzing the remaining equations containing arbitrary functions of one variable we come to 
conclusion that one of them can admit wider invariance groups iff either 

1) it is equivalent to PDE of the form ()6.1() . or 

2) it is equivalent to PDE of the form ()6.15|) . 

Skipping the proof, we present two typical examples. We begin with the equation 

Utt=Uxx + U + G{U;j:)- (6.19) 

This equation is invariant under the four-dimensional algebra {dt, dx, e^du, e~^du) isomorphic to j43.6© 
Ai. Inserting F = u + G{ux) into classifying equation (|2.5|) yields the system of two equations for the 
function G 

h'G' = -h" - 2A, [{h - X)Ux + rx]G' - {h - 2X)G = m - r^x - ^h'u^ - r. 

As we require G" / 0, it follows from the first equation that A = /i' = and the second equation takes 
the form 

{hux + rx)G' - hG = m - r^x - r. 

Upon differentiating the above equation twice with respect to Ux we get {hux + rx)G" = 0. As G" ^ 0, 
relations h = Vx = hold. Hence we conclude that the class of PDEs (|6.19j) does not contain equations 
admitting five-dimensional algebras. 

The system of determining equations for symmetry group of 2'^^'^^^^^^^ equation 

utt = Uxx + w]^'^G{uj), rj = X — kt, u! = r]u~^Ux, k>0, (6.20) 

read as 

(r/"^rw - ri~^rx)Guj - r]~'^rG = ru - Vxx, 

[(A2 - kXi)ri'^uj + ri-^h']G^ - 2(A2 - kXi)r]-^G = -2h'r]-^uj - h" . 

Differentiating the first equation with respect to iv yields 

{7]~^ruj - T]^^rx)G^^ = 0, 

whence in view of inequality G^^^j 7^ we get r = 0. 

Next, differentiating the second equation twice by lo we get 

[(A2 - /A;Ai)7?-3^ + v'^h']G^^^ = 0, 

whence it follows that G^^^^ = 0. Indeed, if this relation does not hold, we have A2 = kXi, h' = and 
operator (|2.4() takes the form 

X{tdt + xdx) + Xi{dt + kdx) + Giudu, A,AiCiGM, A; > 0. 

As the above operator contains at most three arbitrary constants it cannot generate a four-parameter 
Lie transformation group. 

By the above argument we can restrict our considerations to the following class of functions G : 

G = Auj^ + Blu + G, A/ 0,-1, B,GeR. (6.21) 

We can suppose that A ^ —1 in 1)6.21(1 (since otherwise (|6.2U|) belongs to the class of PDEs (|6.1() ). 
Inserting (|6.21j) into the second equation from ((6.211) yields 

2{A + l)r]^h' = B{X2-kXi), if Bh' + ifh" = 2G{X2 - kXi). (6.22) 



If A; > 0, then this system spHts into the following equations (note that h = h{x)): 



h' = 0, B{\2 - kXi) = C{X2 - kXi) = 0. 

Provided \B\ + |C| 7^ 0, there is no way for extending symmetry of equation (|6.2U() . If, on the contrary, 
B = C = 0, then F = Au~^u'^{A ^ 0,-1) and we obtain the equation equivalent to (|6.15|) . Under 
k = system ()6.22() takes the form 

2{A + l)x^h' = X2B, x^Bh' + h"x^ = 2X2C. 

Hence 

h = -]^X2{A+l)-^Bx-' + Ci, CiGM, C=^^^. 

In this case equation (|6.2U|) does admit additional symmetry operator 

dx ^x~^udu 

2(^ + 1) 

but the change of variables 

B 



t = t, X = X, u = \x\ V, 



2{A + 1)' 



reduces it to the form (|6.15jl . 

So equation (|6.2U() admits wider symmetry group iff it is either belongs to the class of H6.1|l or is 
equivalent to (|H.15|) . 

To finalize the procedure of group classification of equations ()2.7|) we need to consider invariant 
equations obtained in the previous section that contain arbitrary functions of two variables. 



6.3.2 Classification of equations with arbitrary functions of two variables. 

In the case under study the standard Lie-Ovsyannikov method is inefficient and we apply our classifi- 
cation algorithm. In order to do this we perform extension of three-dimensional solvable Lie algebras 
to all possible realizations of four-dimensional solvable Lie algebras. The next step will be to check 
which of the obtained realizations are symmetry algebras of nonlinear equations of the form 1)2. 7|) . In 
what follows we use the results of the paper where all inequivalent (within the action of inner 
automorphism group) four-dimensional solvable abstract Lie algebras are given. 

We give full computation details for the case of ^a.e-invariant equations. As shown in the previous 
subsection, there are two inequivalent ^3.6-invariant equations, namely, 

^3.6 = {e~^du, e^du, dt + kd^) 

{k > 0) : F = u + G{r], Ux), rj = x — kt; 

1 3 1 

^3.6 = {\t\'''^du, \t\'^du,tdt + xdx + -udu) : 

3 3 
F = -r'^u + \t\'2G{i,uj), ^ = tx~^, uj = x~^ul. 

According to the algebra ^3.6 is the subalgebra of the following four-dimensional solvable Lie 
algebras: 2A2.2, A3.6 ^1; ^4.2(9 = -1); ^4.8(9 = -|)- 

Algebra 2^2.2- The algebra 2^2.2 = (61,62,63,64) is determined by the following commutation 
relations (note that we give non-zero relations only): 



[ei, 62] = [ei, 64] = [62, 63] = [62, 64] = 0, [61, 62] = 62, [63, 64] = 64. 



It contains a subalgebra ^3.6 = (ei — 63, 62, 64). That is why, we can choose as ei — 63, 62, 64 the basis 
operators of the reahzation of ^3.6- Next, we take as ei + 63 an arbitrary operator of the form ()2.4() 
and require for the commutation relations 

[ei - 63,61 +63] = 0, [61 + 63,62] = 62, [61 + 63,641= 64 (6.23) 

to hold. 

Realization ^3 g. In this case 

ei - 63 = -dt - kdx, 62 = e'^du, 64 = e^du, 
61 + 63 = {\t + Xi)dt + (Ax + A2)5^. + {hu + r)du. 

It follows from dlT^ that 

A = Ai = 0, r = 7 = 7(r/), h = -1. 

Using the change of variables 

t = t, X = X, v = u + A{t]), 
where A = A{rj) is a solutions of equation A2A' + A = —7, we simplify the operator ei + 63 to become 

ei + 63 = adx — udu, a G M. 

Requiring invariance under the above operator yields that a 7^ (otherwise G would be linear in 
Ux)- With this condition we rewrite the invariant equation to become 

G = exp{—a^^ri)H{u)), u> = ex.p{a^^'r])ux- 

Thus we arrive at the following realization of the algebra lAi.2'- 

(6"*5„, 6*5„, dt + kdx, adx - udu) {k>0, a / 0). 

This algebra is admitted by the equation 

utt = Uxx + u + exp(— a~^r/)G(w), rj = x — kt, oj = eiqi{a~^ri)ux- 

If the function G {G^^ 7^ 0) is arbitrary, then the obtained realization is the maximal symmetry 
algebra of the equation under study. What is more, no G exists such that the above equation admits 
a wider invariance algebra. 

Realization ^35- In this case 

ei- e'i = -tdt- xdx-]^udu, 62 = |tr^9„, 

3 

ei + 63 = {\t + \i)dt + {\x + h)dx + (hu + r)5„, 64 = \t\2du. 

It follows from commutation relations 1)6. 23(1 that Ai = A2 = A = 0, /i = —1, r = |i|27(,^), ^ = tx~^. 
Making the change of variables 

t = t, X = X, f = n— |t|27(^) 

we get r = in 61 + 63. Consequently, without loss of generality we can choose 61 + 63 = —udu- 
Requiring for g-invariant equation to admit the operator 61 + 63 yields the equation 2loGu, = —G, 
whence G = \uj\~^H{(^). Consequently, the function F is linear in the variable Ux- This means that 
g does not admit extension to a realization 2A2.2 that can be a symmetry algebra of an equation 
of the form (|T7|) . 



Algebra A^ q © Ai. What we need to do here is to supplement the set of operators ei, 62, 63 forming 
the basis of ^3.6 by the operator 64 of the form ()2.4|) and verify the commutation relations 

[ei, 64] = [62, 64] = [63, 64] = (6.24) 

Realization A'^q. It follows from ()6.24|) that h = X = Xi = 0, r = 7(7/), r] = x — kt in the operator 64 
so that 

64 = aSa; + 7(r/)5„, a e M. 
If Q / 0, then 64 is equivalent to dx- Hence we get two possible realizations of the algebra A3.6 © Ai: 

(e"*5„, e^du,dt,dx); 

{e'^du, e^du,dt + kdx,-f{v)du)- 

Analyzing the above realizations we come to conclusion that the second one cannot be invariance 
algebra of a nonlinear equation of the form (|2.7)) . The first realization is the maximal (if G is an 
arbitrary function) invariance algebra of the equation H6.19jl . 

Realization A^q. It follows from (|6.24|1 that Xi = X2 = h = X = 0, 

so that the operator 64 necessarily takes the form 64 = \t\^^{£,)du- As the straightforward verification 
shows thus obtained realization cannot be invariance algebra of a nonlinear equation of the form H2.7() . 

Algebra ^44.2 {q = —1). We need to supplement the set of operators 61,62,64 forming the basis of 
^3.6 by the operator 63 of the form (|2.4() so that the following commutation relation hold 

[61,63] = [62,63] = 0, [63,64] = 62 + 63. (6.25) 

Realization Ai,.. In this case 



61 = 6 '^du, 62 = 6*9„, 64 = -dt - kdx 

and it follows from 1)6. 25() that coefficients of 63 satisfy equations /i = A = Ai = A2 = 0, the function r 
being a solution of the equation 

T^t + f^T^x = r + 6*. 

This realization cannot be invariance algebra of nonlinear equation of the form (|2.7|1 . 
Realization ^3 g . In this case 



ei = |t| ^du, e2 = \t\2du, = -tdt - xd^ - -udu- 

It follows from (|6.25|) that the coefficients of the operator 63 satisfy equations A = Ai = A2 = /i = 
and the function r is a solution of the equation 

3 , ,3 
trt + xr^ = -r + \t\2. 

Further analysis shows that the so obtained realization cannot be invariance algebra of a nonlinear 
equation of the form (|2.7|) . 

Algebra ^44.8 [q = — ^). We need to supplement the set of operators 61, 63, 64 forming the basis of 
^3.6 by the 62 of the form (|2.4() in order to satisfy the commutation relations 

[61,62] =0, [62,63] = 61, [62,64] =62. (6.26) 



Realization A'^q. In this case 



1 . 1 



ei = e du, 63 = 6 du, 64 = -dt + -kd.. 



X 



and the second commutation relation yields the false equality 1 = 0. 
Realization ^3 g . In this case 

61 = \t\~2du, 62 = \t\^du, 64 = -tdt - xdx - ^udu 

and again the second commutation relation from (|6.26|) cannot be satisfied. 

So that there are no extensions of the realization of ^3 e to a realization of the algebra ^33 (q = 

-\)- 

The remaining equations containing arbitrary functions of two variables are handled in a similar 
way. The results can be summarized as follows 

1) if the functions contained in the equations under study are arbitrary, then the corresponding 
realizations are their maximal invariance algebras, and 

2) except for equation (|6.14l) . all the equations in question do not allow for extension of their 
symmetry. 

Below we give the complete list of invariant equations obtained through group analysis of equations 
with arbitrary functions of two variables. 

6.3.3 Equations invariant under four-dimensional solvable Lie algebras. 

A2.2 © 2Ai-invariant equations 

1) {dx,dt +udu,6^du,6''^du) : F = u + 6^G{io), to = u^^u^; 

2) {^{dt + kd.),6^+'''du,6'^du,d. + udu) {k > 0,7] = X - kt) : 

Ik 

F = {k^ -l)u + 6'^G{u)), uj = e-''{ux-u). 
2^2.2-iiivariant equations 

1) {dt + eudu, dx,6''+^'du, e"-'=*a„) (e = 0, 1; A: > 0) : 
F = (A;2 - l)n + e'*G(w), uj = e~'\ux - u); 

2) {adx - udu, dt + kdx, 6-'du, 6'du) (/c > 0, a > 0) : 

F = u + exp{—Q~^rj)G{uj) , rj = x — kt, to = ex.p{a~^ri)ux. 

A3 3 © ^i-invariant equations 

1) {dt,dx,du,tdu) ■■ F = G{ux). 

Asa © ^i-invariant equations 

1) {du, dx,tdt + xdx + (n + x)du, tdu) : 
F = t~^G{oj), u; = Ua; — In 

2) {dt + udu,dx,tdu,du) : F = e^G{uj), oj = e'^u^; 

3) {x~^du,dx - x~''^{u + \D.\x\)du-,tdt + xdx.,tx'^du) : 
F = 2x~^Ux + + t~"'"x~"'"G(a;), a; = xux + u — In 



yl.3.5 © ^i-invariant equations 



1) (9^, du, tdt + xd^ + udu, tda) : F = t^'^G{ux); 

2) {x~^du, dx - x~^udu,tdt + xd.x,tx^^du) : 

F = -2x~'^u + 2t"^(ua; + x'^u) \n\t{ux + x^^u) 

+t~^{Ux + X~^u)G{lo), UJ = XUx+U. 



^3.6 ® yli-invariant equations 



1) {dx, tdu, tdt + xdx,du) : F = t G{lo), u = t u 

2) {dt,dx,e'du,e-'du) : F = u + G{ux). 

^43.7 © 74i-invariant equations 



1 



1 



1 



1) (exp --(1 - q)t an,exp -(1 - q)t du,dt + -{1 + q)udu,dx) 



{q^O,±l): F=-{1- qfu + exp -(1 + q)t G{u;), 



LO = exp (^-2(1 + 5)*J "a;; 

2) (5,, + + ^(1 + q)udu) 
{q / 0, ±1) : F = - 1)^-2^ + 

a; = |i|5(i-5)«,; 

1 9+1 1+0 1 1 + 9 

3) {\tr~^\(\^du,dx-^x-\du,-q{tdt + xdx),\^\^du) 



{q^O,±l):F 
1+9 



2? 

1 - g2 



ti + 



1 + 



9-1 



g + 1 



2g 



29 G(a;), e = ia;-S a; = l^l'^. xux + 

© ^i-invariant equations 

1) {smtdu,costdu,dt,dx) : F = -u + G{ux). 
^43.9 © 74i-invariant equations 

1) {sm.tdu, cos tdu, dt + qudu, dx){q > 0) : 



-u 



+ e«*G(a;), u = e 



qt 



^4.1 -invariant equations 



1) {du, -tdu, dx, dt - txdu) ■ F = G{u), uj = Ux + -i^; 

2) {du, -tdu, adx + ^t^du, dt + kxdx) (A; > 0, a > 0) : 

F = a-^{x - kt) + G{ux). 



yl.4.2-invariant equations 

1) {\t\^-ydu,\t\''''du,d:„tdt + xd^ + 

/ 0, 1) : F = \qiq - 2)t-^u + \t\\^'i-^^ G{ul 



du) 



Ux - 2\t\2; 



2) {d^, ^/\t\^u, \t\du, tdt + xd^ + + ^^«) 



^4.3-invariant equations 



1) ((9a;, i^i^^u, -|t|2 ln|t|a„,t9t + + ^udu) ■■ 

F = -^t-^u+\t\-^G{cj), Lo = \t\^Ux; 

2) {dx,tdu, du, tdt + xdx) : F = t~'^G{u), uj = tu^; 

3) {e'^'du, dt + kudu, Pdx + te^'du, e'^'du) (fc 7^ 0, ^ > 0) : 
F = k^u + 2A;/3-^xe** + e^*G{u), u = 6"*^*^^; 

4) {e'+^'du, e^'du, a{dx + n5„) + 2He^5„, -^{dt + kd^)) 

(a 7^ 0, A; > 0) : F = {k^ - \)u - Ak'^a'^rie'^ + e''G(a;), 
LO = e~^{ux — u),rj = X — kt. 



A4.4-invariant equations 



1) {\t\^du, -\t\nn \t\du, dx, tdt + xdx + 
1 



-u - x\t\2 ln\t\ 



du): 



F = ^t-^u + \t\-^G{Lj), u = \t\-'^ux + ^ In^ \t\. 



A4,5-invariant equations 



1) {dx, \tr-''du, \t\'-'^+''du,tdt + xdx + mudu) 

(m/ -(l + a),-+Q;;a/0) : 

F = {m- a){m -a- l)t-\ + \tr-^Giuj), u = 



^4.6-invariant equations 



1) {dx, \t\^sui{q-^hi\t\)du, \t\^ cos{q-^hi\t\)du,qtdt + qxd^ 



F = - Q + + |t|«"'(f-i9)G(a;), u = \t\i~'ih-p) 

^4.7-invariant equations 

1) {du, -tdu, dt + kdx,tdt + xdx + (2u - h''^ du) {k>0): 
F = — In + G(ti;), oj = r]~^Ux, r] = x — kt. 



^4.8"iiivariant equations 



1) {dt + eudu, d,, e'duMdu) (e = 0; 1) : 
F = -u + e'*G(w), uj = e~^\ux - u)- 

2) (Ixp-'^a^, du t\xr-'^du,tdt + xd^ + mudu) (g / O, m G M) : 
F = -{m -q){m- q - l)x"^n + |x|""^G(a;), 



3) (9t + kdx, du,tdu, tdt + xd^ + qudu) {k > 0,q eR) : 
F = \t]\'}-^G{u;), lo = M^-'^u^, i] = x-kt; 

4) {x~^du, dt + dx- x'^udu, tx'^du, tdt + xd^) : 

F = 2x~^Ux + x~^{t — x)^^G{uj), oj = xux + u; 

5) {du, -tdu, dt + kdx + udu, ad^ + udu) (a / 0, A; > 0) : 

F = exp(a^^r/ + t)G{uj), oj = exp(— a~"'^r/ — i] = x — kt. 

^4.10-iiivariant equations 

1) (sin tdu , cos tdu , dx + udu , dt + kdx) {k > 0) : 
F = —u + e^G{uj), ijj = e~^Ux, rj = x — kt. 

In the above formulas G = G{uj) is an arbitrary function satisfying the condition -Fu^u^ 7^ 0. 

7 Symmetry reduction and solutions of nonlinear wave equations 

Among various apphcations of Lie symmetry groups the most prominent and remarkable one is a 
possibility to construct exact solutions of nonlinear PDEs. The basic idea is reducing multi-dimensional 
differential equations to ordinary differential equations via special ansatzes (invariant solutions). A 
regular (but not the only) way to derive those ansatzes is to utilize symmetry group admitted by the 
equation under study (for more details see, e.g., ^lEl)- Though the obtained ordinary differential 
equations are, as a rule, nonlinear, they possess in many cases a residual symmetry allowing for 
constructing their general or particular solutions. Inserting the latter into the corresponding ansatz 
yields the exact solution of initial nonlinear multi-dimensional PDE. This method is often referred to 
in the literature as symmetry reduction of PDEs. 

We apply the symmetry reduction approach to derive the families of exact solutions of nonlinear 
wave equations ()2.7() having the richest symmetry properties. 

To perform reduction of PDEs (|2.7|) to ordinary differential equations we need to obtain all in- 
equivalent one-dimensional subalgebras of the symmetry algebras of the equations under study. What 
is more, basis operators of the said one-dimensional algebras 

T(t, X, u)dt + i{t, X, u)dx + r]{t, x, u)du, 

have to obey the following restriction jl4j : 

|T| + |e|^0 (7.1) 

in some open domain Q of the space y = x of independent = {t, x) and dependent M} = {u) 
variables. 



As we proved in the previous sections, equations 



Utt = Uxx - u~^ul; (7.2) 

Utt = ^^xx + e""; (7.3) 

Utt = Uxx + mlnluxl {m > 0); (7.4) 

Utt = Uxx + l^^xl'' {k + 0,1). (7.5) 



enjoy the highest symmetry properties amongst PDEs of the form ()2.7() . 

The first step of symmetry reduction algorithm is classifying one-dimensional subalgebras of the 
invariance algebras of the above equations taking into account constraint 1)7.11) . To classify one- 
dimensional subalgebras we utilize the method suggested in ||45j and the lists of one-dimensional 
subalgebras of four-dimensional subalgebras obtained in 0Sj . 

Equation (|7.2|) admits the algebra 

A\ = (.43.3 ©^i)^) (es), 

where ^3.3 = (61,62,63) = {udu,dx,xudu) , Ai = (64) = (9i),e5 = tdt + xdx- 

In what follows we will need the commutation relations of the basis operators of the algebra 
j43.3 © Ai with the operator 65: 

[61,65] =0, [62,65] =62, [63,65] = -63, [64,65] =64. 

According to one-dimensional subalgebras of ^3.3 © Ai defined within action of inner automor- 
phism group of this algebra read as (ei), (61-1-064), (64), (62), (62 -|- 064), (63), (63 -F 064), {e2 + Pe3), 
(62 + /363 -|- 064) {a, (3 / 0). The above subalgebras can be further simplified by using transformation 
group generated by the operator 65. For example, using the Campbell-Hausdorff formula we transform 
61 -|- 064 as follows: 

exp(065)(6i -|- 064) exp(— ^65) = 61 -|- 06^64. 

Consequently, putting 6 = — ln|a| we simplify 61 -|- 064 to become ei it 64. Similarly, we prove that 
we can put a = ±1 in 63 -|- a64 and /3 = ±1 in 62 -|- Pe^, 63 -|- /?63 -|- a64. 

To complete classification of one-dimensional subalgebras we have to describe all inequivalent 
subalgebras with non-zero projection on the basis operator 65, i.e., subalgebras of the form 

A = 65 -I- ai6i -I- 0262 -I- 0363 + 0464, ai, a2, "3, 04 G R. (7.6) 

Utilizing the automorphism exp(0e4) with properly chosen 9 we have 04 = in ()7.H|) . Next, applying 
transformation exp(0i62 -|- ^263) to operator 1)7. 6|) reduces it to one of the following operators: 65, 65 -|- 
a6i (a / 0). 

So the list of one-dimensional subalgebras of the five-dimensional algebra A^ determined up to the 
action of inner automorphism group is exhausted by the following algebras: (61), (61 it 64), (64), (62), 
(62 -|- 064), (63), (63 ± 64), (62 ± 63), (62 ±63-1- 064), (65), (65 -|- aei) (q / 0). By direct verification we 
prove that the basis operators of the algebras (61), (63) do not satisfy condition ()7.1j) . 

Finally, we make use of the fact that the discrete groups of transformations 

t = —t, X = X, V = u; 
t = t, X = —X, V = u; 
t = —t, X = —X, V = u, 

also belong to the equivalence group of (|7.2)) . Using the above transformations enables to further 
simplify the optimal system of inequivalent subalgebras 



(61-^64), (64), (62), (62-^064), (63-^64), 

(62 ± 63), (62 ±63 + 064), (65), (65 + aei) (a > 0). 



(7.7) 



The second step of the method of symmetry reduction is constructing the complete set of invariants 
f{t, X, u) for each inequivalent one-dimensional subalgebra. As a typical example, we consider the case 
of subalgebra (ei + 64). To construct its invariants we need to integrate the first-order PDE 

(ei -|- 64) o F{t, x,u) = 



or 

uFu + Ft = 0. 

The complete set of first integrals of the above equation reads as wi = x, uj2 = e~*n. Hence we get 
the general form of invariant solution (ansatz) uj2 = fii-oi)- Solving this equation with respect to u we 
finally have 

u = e^ip{x). (7.8) 
Inserting (|7.8|) into (|7.2|) yields ordinary differential equation for unknown function (p 



which is easily integrated 



if = exp 



+ Cix + C2 



Inserting the so obtained expression for (j) into ansatz (|7.8() yields the explicit form of invariant solution 
of equation ()7.2|) 



u = exp 



t + -x^ + Cix + C2 



, Ci, C2 € M. 



Analysis of the remaining subalgebras from (|7.7j) yields the following results. Invariant solutions 
of H7.2() that correspond to the algebras (64), (62), (e2 + 064) (q > 0) read as 



u 
u 
u 



exp[Cix -I- C2], 
Cit + C2, 

Calt-ax + Cil^-"', 



where Ci,C2 G M. 

The ansatz corresponding to the algebra (63 -|- 64) has the form 

u = e*'^<p{x). 

It reduces (|7.2j) to ordinary differential equation 

if — (y9 [if ) — X(/9 = L), 

whose general solution has the form 

" 1 



Lp = exp 



12 



-x" + Cix + C2 



, Ci, C2 G M. 



Inserting the above expression into the corresponding ansatz we get another family of exact solutions 
of nonlinear wave equation (|7.2)) 



u = exp 



tx + ^x'^ + Cix + C2 



Ci,C2 G M. 



A similar analysis applied to the algebra (e2 + 663) (e = ±1) yields two families of exact solutions 
(El: 

u = Ci exp ^— -x^^ cos(t + C2), e = — 1; 
u = Ci exp ^2^^^ cosh(t + C2), e = 1; 

where Ci,C2 G M. 

The ansatz invariant under the algebra (e2 + ees + 064) (e = ±1; a > 0) reads as 

u = exp ^2^^^^ fiv)j rj = t — ax. 

It reduces PDE H7.2|l to equation 

(a^ - l)ip" - a^ip~\^'f + = 0, 
whose general solution is given by one of formulas below 
u = exp(Ci±?7), a = l,e = l; 



U 



C2 
C2 



COS I Ci + 



7] 



cosh Ci + 



a2-l 
V 



1 - a' 



a > 0, a ^ 1, e = — 1; 



a > 0, a / 1, e = 1; 



where Ci, C2 G M. The corresponding invariant solutions of (|7.2jl have the form 



u 



u 



u 



exp 



Ci + -x±{t- x) 



C2exp ( --x^ 



C2exp ( --x^ 



cos I Ci + 



cosh Ci + 



t — ax 
t — ax 



1-a' 



1 



a > 0, a ^ 1; 



, a > 0, a / 1; 



where Ci,C2 G M. 

The algebras (es), (es + aei) (a > 0) give rise to the so called auto-model (J^) solutions of H7.2() . 
Solution invariant under (es) reads as 



exp 



|/ (ci(r/2-l) + i(r?2-l)ln 



1 +?? 



1 — ?7 



2^ 



c^?? + C2 ^ , 



C = te-\ Ci,C2GM. 
Solution invariant under the algebra (es + aei) {a > 0) has the form 

n=|t|"(^(^), ^ = te~^ 
Inserting this expression into H7.2() yields ordinary differential equation 



If a = 1, then the general solution of the above equation is of the form 



Cirf + -rj In 



1 + r? 



1 — 7/ 



-1 



dr] + Ci 



if = exp 

i = tx'^, Ci,C2GM. 
Provided a > 0, a / 1, the change of variables 

ip = expf, f = fiO 

reduces the equation for (p to become 

fif - 1)/" - fifT + mf - o^)f - c^io^ - 1) = 0. 

Now making the change of variables, f = g, g = g{(,) we arrive at the Riccatti equation 

fif - 1)5' = eV - 2e(c' - a)g + «(a - 1). 



(7.9) 



In view of the above we have the two families of exact solutions of nonlinear wave equation 1)7.2(1 : 



u 



u 



t exp 



C1I] +-r?ln 



1 + ry 



1 — ?7 



- 7] 



dr] + C2 



Itl^exp [ g{r])dr] + Ci 



where a > 0, a / 1, ^ = tx ^ ,Ci & E., and the function g = g{^) is a solution of (|7.9j) . 

Now we turn to nonlinear wave equation (|7.3|) . Its maximal symmetry algebra is A"^ = (^3.3©^!)^ 
(es), where ^3.3 = (61,62,63) = {du,dt,tdu), A\ = (64) = {dx), 65 = tdt + xdx + {u - x)du- What is 
more, the basis operators of the algebra ^43.3 © Ai obey the following commutation relations with 65: 

[61,65] =61, [62,65] =62, [63,65] =0, [64,65] = 64 - 61. 

Using these and taking into account condition (|7.1|) we obtain the following system of optimal one- 
dimensional subalgebras of the algebra A^: 

(64), (62), (62 + 064), (63 + 064), (62 ± 63), 

(62 + /363 + 064), (65), (65 + Q63) (a > 0, /? 7^ 0). 

Skipping the intermediate computations we give the final result, the families of exact solutions of 
nonlinear equation (|7.3|1 : 

(64) : n = Cit + C2, Ci,C72GK; 

(62) :n = (x + Ci)[l-ln|x + C7i|] + C2, ^1,(^2 GM; 

(62 + 064) : u = {x — at) In |1 — a^l 

-(x - at + Ci)[ln \x- at + Ci\ - 1] + C2, 

a>0,a/l, Ci,C72GM; 
(63 + 064) : u = a~^tx + o? exp(a~^t) + C\t + C2, a > 0, Ci, C2 G M; 



(62 + ees) : u = ^et^ + ip{x), e = ±1, ip' = y{x), 
y-ln\e-ey\ = ex + Ci, Ci G M; 

{62 + ^63 + 64} ■.u = ^l3t^ + ix-t)ln\p\+Ci, CieR; 

{€2 + Pes + ae^) : u = ^Pt'^ + ip{r]), rj = x - at, (f = y'{r]), 

y-ln|/3-e^| =^(1-02)^^ + ^1, a > 0, a / 1, /3 / 0, Ci G M; 
(es) : n = -xln|x| +Cit + x, Ci G M; 
(es + QCs) : u = (at — x) In |x| + xf{^), ^ = tx^^, a > 0, 
- l)ip" + exp[-^(^' + ip + aC-l] = l + aC. 

Equation (|7.4j) admits the algebra = (A3. 5 © j4i)^ (es), and A3. 5 = (61,62,63) = {du,dt,tdu), 
Ai = (64) = (dx), 65 = t9f + xdx + (2ti + ^nit'^) du- Its optimal system of one-dimensional subalgebras 
has the form 

(ei + 64), (64), (62), (62 + 064), (63 + 064), (65) (q > 0). 

Below we give the list of exact solutions of (|7.4j) invariant under the above subalgebras. 

(61 + 64) : u = X + Cit + C2, Ci, C2 G M; 
(64) :u = Cit + C2,Ci,C2£M; 

(62) : u = Lp{x), If' = /(x), / - — |— 7 = -mx + Ci, Ci G M; 



In I/I 

(62 + 064) : n = X — t + Ci, Ci G M, if a = 1; 
u = ip{r]), r] = x-at, ip' = /(r/), 

/df Tfi 
= n{x -at) + Ci, Ci G M, if a > 0, a / 1; 
In I/I l-a^ 

(63 + 064) : u = a^^tx + —mt^{l + Ina) + 2"^^^ ^^'^ 1^1 ~ 2) 
+Cit + C2,a> 0,Ci,C2 G M. 

(65) : li = -w? In |t| + 4^99, where (p = p{^) = tx^^) satisfy equation 
^2(^2 _ 1)^// + 2e(^2 - 2)(^' - + mln |eV'l - = 0. 

Note that the particular solution of reduced equation (/? = exp(|) gives rise to the following 
invariant solution of (|7.4j) : 

u = 2"^^^ ^n |t| — tx exp(— ). 

Equation 1)7.5^ is invariant under the algebra A| = (^3.5 © Ai)l^ (65), and ^3.5 = (61,62,63) = 
{du,dt,tdu), A\ = (64) = {dx), 65 = tdt + xdx + ^E^udu, fc / 0; 1. The optimal system of one- 
dimensional subalgebras of this algebra reads as (65)5(61 + 64), (64), (62), (62 + 064), (63 + 64)5(62 + 63), 
(62 + 63 + 064) (a > 0). If /c = 2, then the above list should include also the algebra (65 + 061) (a > 0). 



Below we give exact solutions of nonlinear wave equation (|7.5jl invariant under the above algebras. 

(ei +64) : n = + x + Cit + C2, Ci,C2eM; 
{ei) : u = Cit + C2, C7i,C2GM; 

(62) : -u = (2 - k)-^{Ci + {k- l)x)^^ + C2, Ci, C2 G M, 
if A: / 2; 

n = C2ln|x - Cil, Ci,C2eM, if /c = 2; 
(e2 + 064) : u = C, C G M, if a = 1; 

u = {l-a^)\n\Ci-x + at\+C2, Ci,C2eR, 
if k = 2,a > 0,a ^ I; 

1-k / l-k\^ 



u 



|x-Qt + Cl|l-'= +C2, Cl,C2 G 



2 - V"^ - 1, 
if / 0,l,2;a > 0,a / 1; 

(es + 64) : u = + t(ln |t| - 1) + Cit + C72, Ci, C2 G M, 
if /c = -1; 

u = - In |t| + Cit + C2, Ci, C2 G M, if A; = -2; 
n = tx+(A;2 + 3A; + 2)^i|t|'^+2^C7it + C72, ^1,^2 G 

if /c / 0,1,-1,-2; 

1 /"^ 

(e2 + ees) : u = ^^f^V, f = fix) = / /(??)dr/ + C2, 

where / is defined by 
df 



e-l/l' 



7? + C7i,Ci,C2 GM, e = ±l; 



1 

(62 + eea + ae4) : u = -et^ + if, ip = (p{ri) = / f{z)dz + C2, 



r] = X — at, where / is defined by 



e-l/l' 



(1 - a^)-^z + Ci, Ci, C2 G M, e = ±1, a > 0; 



fc-2 _, 

(es) ■.u = \t\ '="^93(0, = tx~ , where is defined by 



k-2 
k-l 



+ T^^^ = 0, if A;/ 0,1, 2; 



Ci(l-r?2) + -(1-^2)1^ 



1 + r? 



1 — T] 



dr] + C2, 



Ci,C2 G = if A; = 2; 

(es + aei) : u = a In |t| + ip = (^(^) = y" /(r/)(ir/ + C, 

^ = C G M, / = /(?]) is a solution of Riccatti equation 
rj^irj"^ - 1)/' + + 2r]^f + a = 0, a > 0, if A; = 2. 



8 Concluding remarks 



Let us briefly summarize the results obtained in this paper. 

We prove that the problem of group classification of the general quasi-linear hyperbolic type 
equation reduces to classifying equations of more specific forms 

I. Utt = Uxx + F{t, 

II. Utt = Uxx + git,x,u)ux + f{t,x,u), 5« / 0; 

III. Utx = g{t,x)Ux + f{t,X,u), 9x7^0, fuu^O; 

IV. utx = fit,x,u), fuu^O- 

The cases of PDEs that are essentially nonlinear in Ux (the class of PDEs I) and either linear in Ux or 
do not depend on Ux (the classes II - IV) need to be considered separately. 

If we denote as T>£ the set of PDEs II - III, then the results of application of our algorithm for 
group classification of equations I - IV can be summarized as follows. 

1) We perform complete group classification of the class T>£. We prove that the Liouville equation 
has the highest symmetry properties among equations from DS. Next, we prove that the only 
equation belonging to this class and admitting the four-dimensional invariance algebra is the 
nonlinear d'Alembert equations. It is established that there are twelve inequivalent equations 
from T>£ invariant under three-dimensional Lie algebras. We give the lists of all inequivalent 
equations from DS that admit one- and two-dimensional symmetry algebras. 

2) We have studied the structure of invariance algebras admitted by nonlinear equations from the 
class I. It is proved, in particular, that their invariance algebras are necessarily solvable. 

3) We perform complete group classification of nonlinear equations from the class of PDEs I. We 
prove that the highest symmetry algebras admitted by those equations are five-dimensional 
and construct all inequivalent classes of equations invariant with respect to five-dimensional Lie 
algebras. We also construct all inequivalent equations of the form I admitting one-, two-, three- 
and four-dimensional Lie algebras. 

The results of group classification of the class of nonlinear wave equations are utilized for 

constructing their explicit solutions. Namely, we perform symmetry reduction of all equations 
admitting five-dimensional symmetry algebras to ordinary differential equations and constructed multi- 
parameter families of their exact solutions. 
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